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ABSTRACT 

Integral field spectrographs are major instruments to study the mechanisms involved in 
the formation and the evolution of early galaxies. When combined with multi-object 
spectroscopy, those spectrographs can behave as machines used to derive physical 
parameters of galaxies during their formation process. Up to now, there is only one 
available spectrograph with multiple integral field units, e.g. FLAMES/GIRAFFE 
on the VLT. However, current ground based instruments suffer from a degradation 
of their spatial resolution due to atmospheric turbulence. In this article we describe 
the performance of FALCON, an original concept of a new generation multi-object 
integral field spectrograph with adaptive optics for the ESO Very Large Telescope. 
The goal of FALCON is to combine high angular resolution (0.25 arcsec) and high 
spectral resolution {R > 5000) in J and H bands over a wide field of view (10 x 
10 arcmin^) in the VLT Nasmyth focal plane. However, instead of correcting the 
whole field, FALCON will use multi-object adaptive optics (MOAO) to perform locally 
on each scientific target the adaptive optics correction. This requires then to use 
atmospheric tomography in order to use suitable natural guide stars for wavefront 
sensing. We will show that merging MOAO and atmospheric tomography allows us to 
determine the internal kinematics of distant galaxies up to z « 2 with a sky coverage 
of 50%, even for objects observed near the galactic pole. The application of such a 
concept to Extremely Large Telescopes seems therefore to be a very promising way to 
study galaxy evolution from z = 1 to redshifts as high as z = 7. 

Key words: galaxies: high redshift, galaxies: kinematics and dynamics, instrumen- 
tation: adaptive optics, instrumentation: spectrographs, methods: numerical 



1 INTRODUCTION 

Thanks to the Hubble Space Telescope, astronomers have 
been able to determine the morphology of distant galax- 
ies located at z > 0.5, showing that galaxies in the past 
were mos tly irregular and smaller than those in the local 
unive rse l|Abraham fc van den Berghl 200 ll: van den Berghl 
2002|'l. and had hi gher merging rates i Le Fevre et al.l 2000l : 
Bundv et al1 l2004'). The use of spatially resolved color maps 
has also allowed to highlight regions of star formation as 
well as showing that galaxies cores were bluer than toda y's 
galaxies bulges (|Abraham et al.lll999l : IZheng et al.ll2004r ). 
Spectroscopic studies have also taught us that star for- 
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mation rates (SFR) were higher in distant galaxies, in par- 
ticular in Luminous Infrared Galaxies (LIRGs) where SFR 
can reach 100 Mq/ut, and that the density of star forma- 
tion was much higher a t z = 1 than today (|Madau et al.l 
1 19961 : iFlores et al.ril999l ) . Therefore, when the re-emission 
from the dust in the IR is taken into account, a simple inte- 
gration of the global star formation history shows that half 
of the present-day stars have been forming s ince z ~ 1 — 1.5 
(|Dickinson et al.ll2003l : iHammer et~aLll2005l ). 

However, the dynamical information of distant galaxies, 
crucial for their studies, still is not well known. Indeed, 
spatially resolved color maps allow us to show where the 
star formation occurs, but not to know how masses and 
gas are distributed in those galaxies, and neither to know 
the physical and chemical properties of the gas. Also, the 
question of the evolution of the fraction of barred galaxies is 
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still c ontroversial today llAbraham et al] Il999l : ISheth et al.l 
I2OO3I : IZheng et all 120051 ). At last, the way galaxies in 
merging systems were assembling their masses in the past 
is still mysterious, as well as the influence of dark matter 
into those exchanges. Kinematics and chemistry of galaxies 
up to 2: « 2 is therefore required in order to evaluate the 
velocity fields of their main components (bulges and disks), 
to determine how important are the respective roles of the 
merging phenomenon and the dark matter, and finally to 
establish the physical origin of the Hubble Sequence. At 
these redshifts, the important features such as the Ha 
emission lines or the stellar absorption lines are redshifted 
in the near-infrared, in J and H bands. This means that 
integral field spectroscopy (IFS) is required, with a high 
spectral resolution (3000 ^ i? 15000), but also with 
a good spatial sampling (1-2 kpc) of the velocity field, 
implying a spatial resolution (0.15 — 0.25 arcaec) better 
than the atmospheric seeing. Moreover, there are others 
important issues such as the the field of view (FoV) and 
the multiplex capability. Indeed, it is important to observe 
galaxies on scales greater than the correlation length (4 to 
9 Mpc), leading to a minimal FoV of 100 square arcmin. 
Observations of several cosmological fields will then allow 
us not to be sensitive to cosmic variance effects, and 
a multiplex capability is then required so that several 
objects can be observed simultaneously, providing a gain in 
exposure time. 

In order to reach the req uired image q uality, the solu- 
tion is to use adaptive optics (|Roddierlll999l ). Adaptive Op- 
tics (AO) allow us to compensate in real-time for the wave- 
front's distorsion caused by the turbulence in the Earth's 
atmosphere, providing to the telescopes their nominal an- 
gular resolution (diffraction-limited imaging). However, the 
FoV on which this resolution is achieved is very low. In- 
deed, all the AO systems built before now work with the 
classical AO method, where a guide star (natural or artifi- 
cial), bright enough and close enough to the science target, is 
used to properly measure the incoming turbulent wavefront 
and compensate for it. In these conditions, the compensated 
FoV around the guide star is equal to a few times the iso- 
planatic patch 9o, eg typically 10 to 20 arcsec in H band. 
When applied to extragalactic astronomy, this makes the 
use of natural guide stars impossible for AO, at least under 
its "classical" form, as the sky coverage would then be lower 
than 3% due to the lack of sufficiently bright guide stars, 
even at a galactic latitude of 30° . 

Several concepts have therefore been proposed in the 
last years to increase the sky coverage a nd the corrected 
FoV of AO systems. Laser Guide Star^ (|Fov fc Labevrid 
1 19851 ) create an artificial star in the sodium layer of 
the atmosphere at an altitude of 90 kms thanks to a 
laser. However this method suffers f rom the cone effect 
and the tilt determinat ion problem l|Tallon fc Fovl Il990l : 
iRigaut fc GendronI \l99^ . meaning that a natural GS is 
always required close to the science target to measure the 
low orders of the turbulent wavefront (tip-tilts). Moreover, 
in the case of multi-objects instruments such as those re- 
quired for extragalactic astronomy, one LGS is required per 



scientific target as well as a dedicated wavefront sensor to 
measure the low orders. This can dramatically increase the 
cost of such a multi-objec t inst r ument. Multi-C onjugated 
Adaptive Opticfl (,Dic5 1 19751: j Johnston fc~Wels h 199i; 



IFusco et al.ll200ll : iLe Louarnll2002l ) have also been proposed 
to improve the FoV of AO systems, using several deformable 
mirrors conjugated to the dominant turbulent layers. But 
the compensated FoV provided by such systems remains 
generally insufficient for extragalactic studies. As an exam- 
ple, the compensated FoV provided by the Gemini-South 
MCAO system or the MAD syst em for the VLT will 
have a diameter of o nly 2 arcmin IjEUerbroek et al. I I2OO3I : 
[Marchctti ct al.' 200>|). At last. Ground Layers Adaptive 
OpticfEl (.Rigaut 20021 ) have also been proposed to widen 
the FoV of AO systems. But this approach assumes that 
most of the turbulence is located close to the telescope 
pupil. Moreover, the corrected FoV remains generally 
insufficient for extragalactic studies, and most of the GLAO 
systems studied until today al ways use at least on e LGS 
to reach decent performance l|Le Louarn fc HubinI |2004| : 
iMorris et aLlbOO^ ). 

Such considerations initiat ed therefore the FALCON 
concept iHammer et all 120021 '). FALCON {Fiber optics 
spectrograph with Adaptive optics on Large fields to Correct 
at Optical and Near-infrared) is a concept of a new gener- 
ation multi-object integral field unit (IFU) spectrograph, 
working with adaptive optics over a very wide field of view 
(10 X 10 arcmin^) at the Nasmyth focus of the VLT and 
suitable for extragalactic studies. But instead of correcting 
the whole FoV, the AO correction is only performed on 
the regions of interest, eg the spectroscopic IFUs observing 
galaxies, thanks to a Deformable Mirror (DM) conjugated 
to the pupil for each IFU. However, with such an architec- 
ture, the sky coverage problem of classical AO remains. We 
have focused on this problem, which is solved by using a 
new type of wavefront reconstruction, inspired from MCAO 
techniques (|Fusco et al.ll200l] : iTokovinin et al.ll200ll ). where 
the WFS measurements from several off-axis NGS around 
each galaxy are used to extrapolate the on-axis galaxy 
wavefront in the pupil. Preli minary perform a nce o f this 
concept have been given in lAssemat et al.l l|2004 ) and 
iHaminer et"aLl l|2004 ). 

In this paper, we will show more accurately the expected 
performance achieved by FALCON thanks to the combina- 
tion of 3D spectroscopy with AO. Firstly, we will present 
the technical specifications of FALCON, derived from the 
higher level scientific requirements. We will then show the 
gain AO can bring to Integral Field Spectroscopy. Following 
that, we will accurately describe the principle of the FAL- 
CON AO system. The next section will show the expected 
performance of such a system in the case of median atmo- 
spheric conditions for the Cerro Paranal. The discussion will 
then compare these performance with those from a GLAO 
system working in the same conditions, and will deal with 
additional source of errors whose influence remains to be 
quantified. 
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2 HIGH LEVEL SPECIFICATIONS 
2.1 Wavelength range 

As explained in the introduction, the dynamical informa- 
tion is required to improve our knowledge of the mech- 
anisms responsible for the formation of galaxies. The ve- 
locity field of the galaxies has therefore to be probed by 
measuring the redshifts of emission lines, that we choose 
to peak in the 1.00 — 1.85 fim wavelength range (cover- 
ing J and H bands), thus avoiding the thermal domain 
(A > 1.95 microns) where the instrumental thermal back- 
ground will dominate the noise. Several emission lines can 
then be used, such as [Oil] (A = 3727A), Hf3 (A = 4861 A), 
[Oil I] (A = 5007A) and in particular Ha (A = 6563A). In- 
deed, this latte r suffers less from extinction than the other 
emission lines l|Liang et al.ll2004[ ). and can be used to map 
the dynamical information up to z = 1.8. The use of shorter 
wavelength emission lines will allow to observe galaxies up 
to g = 2.5. howe ver extinction may then be a problem 
(|Liang et al.ll2004l l. 



2.2 Angular resolution 

Morphological studies of galaxies located into the Hubble 
Deep Fields and the GOODS fields have shown that galaxies 
with z ^ 1 had half-light radii ri/2 smaller th an 0.5 arcsec 
l|Marleau fc Simard|[l998l : iFerguson et~ aL''2Q04'') , with anav- 
erage (fi/2) = 0.25 arcsec at z = 2 ( jBouwens et al.l [20041 ). 
As we want to be study the dynamical processes occuring 
within the galaxies, we must be able to resolve their half- 
light radius, up to z = 2. Therefore, an angular resolution 
of at least 0.25 arcsec is required, implying a pixel sampling 
of 0.125 arcsec. Such a spatial resolution is far beyond abil- 
ities of current ground-based, seeing-limited, integral fields 
spectrographs, as it is definitely better than the atmospheric 
seeing. This strongly suggests the use of adaptive optics, as 
we will see later in section [3] 



2.3 Spectral resolution 

We have to observe in the near infra-red (NIR) to fulfill 
our scientific requirements. However, the spectral window we 
want to explore is cro wded with thin , intense terrestrial OH 
emission lines ( i Maihara et al.l 1 19931 : iRousselot et al. I I2OOOI : 
iHanuschikl bOoiT Observations have therefore to be per- 
formed between these lines, implying to separate them. This 
leads therefore to a minimum spectral resolution R — 5000. 
Considering now the impact on the determination of the 
galaxy dynamics, this is equivalent to resolve velocity dis- 
persions a ^ 25.5 km/s a.t z — 0. 

However, for such an instrument dedicated to the ob- 
servation of distant galaxies, the optimal spectral resolution 
is a compromise between the scientific goals on one hand, 
and the drastic loss in spectral signal-to-noise ratio with 
increasing R on the other. As an example, spectral resolu- 
tions 7? ^ 15000 will allow us to resolve velocity dispersions 
a > 8/(1 + 2) km/s, i.e. greater than 4 km/s for a galaxy 
located at z = 1. A value of 15000 seems therefore to be a 
reasonable upper limit to consider for the spectral resolution 
R. Recent observa tions of z = 0.6 galaxies with GIRAFFE 
jFlores et al.ll2006l ) have in fact shown that a slightly lower 



spectral resolution R = 10000 is definitely required to probe 
precisely the kinematics of distant galaxies with various mor- 
phologies, including compact galaxies. Moreover, a spectral 
resolution R — 10000 also allows us to resolve the [0//J3727 
doublet, and then to retrieve the kinematics of galaxies with 
higher redshift. 



2.4 Field of view and multiplex capability 

Observations of Lyman Break Galaxies with z ^ 2 show 
a strong spat ial clustering, w ith correlation lengths greater 
than 4 Mpc (|Giavaliscdl2002l ). Similar results are also con- 
firmed by cosmic web simulation codes, based on A.CDM 
theory ( Hatton et al.l |2003| ). The observation of distant 
galaxies over scales greater than this correlation length 
should therefore be achievable, in order to avoid statisti- 
cal biases due to the cosmic variance effects. As an example, 
such ef fects have been obs erved on the WFPC2 Hubble Deep 
Fields iLabbe et al.ll2003l ). where the FoV (^ 6 arcmin"^) is 
unsufficient. It is then required to encompass the correla- 
tion lengths for all redshifts, leading to a instrument with a 
minimum FoV of 100 arcmin^ . 

The instrument's FoV is also linked to the number of 
IFU s as well as t he w ished multiplex gain. As an exam- 
ple, ISteidel et al.l (|2004h found a density of 3.8 galaxies per 
arcmin^, for objects with redshifts 1.5 < z < 2.0 and mag- 
nitudes R ^ 25.5. If we assume now that we only observe 
objects with magnitudes R 23.5 or even R ^ 23.0 (be- 
cause of SNR issues), we find a lower source density, com- 
prised between 0.31 and 0.84 galaxies per arcmin'^. From 
our past experience o n ESO large program with GIRAFFE 
(jHammer et al.lll999l '). we think that 15 to 30 IFUs should 
be a reasonable number while providing an efficient gain 
in terms of exposure time. However, in order to optimise 
the exposure time, it is required to observe objects in the 
same magnitude range, and the appropriate density to con- 
sider becomes the number of objects per magnitude and per 
arcmin'^. Moreover, a non negligible fraction of the objects 
will see their emission lines match with atmospheric OH 
lines, making them impossible to observe. As a consequence, 
this might make the source density definitely smaller than 
the values written above, and low enough to warrant a wide 
FoV, equal to 10 x 10 arcmin^ . The VLT Nasmyth focal 
plane is therefore ideally suited for a field of view, as the 
vignetting on a 10 arcmin diameter FoV is less than 1%. 
Spreading several IFUs over the Nasmyth focus will there- 
fore allow us to perform the simultaneous 3D spectroscopy 
of several galaxies with similar redshifts, thus providing a 
very important gain in exposure time. 



2.5 Conclusion 

In this section we have given high level specifications for the 
FALCON instrument, derived from scientific requirements. 
As it has been shown, a multi-object spectrograpy with both 
high spatial {FWHM ^ 0.25 arcsec) and spectral {R ^ 
5000) resolutions is required, moreover over a wide field of 
view (FoV^ 100 arcmin^). 

We would especially like to insist now on the spatial 
resolution we require (0.25 arcsec), which is not achievable 
with current seeing limited instruments such as GIRAFFE. 
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This requires the image quality to be enhanced, and the best 
way to achieve this, in particular in terms of signal to noise 
ratio (SNR), is to use Adaptive Optics. We therefore show 
in the following section the gain brought by Adaptive Optics 
to Integral Field Spectroscopy. 



3 3D SPECTROSCOPY AND ADAPTIVE 
OPTICS 

3.1 Introduction: why use AO with 3D 
Spectroscopy 

Adaptive optics (AO) is a technique that allows us to re- 
store in real-time the flatness of a wavefront distorded by 
the effects of the atmopheric turbulence, thus improving the 
image quality at the focal plane of an instrument. For classi- 
cal imaging applications, the main interest of AO is to allow 
us the partial restoration of spatial frequencies in the images 
up to the telescope's cut-off frequency (D/X), far beyond the 
turbulent cut-off frequency. This is the consequence of the 
morphology of compensated point-spread function (PSF), 
which is quite complex: it cannot only be summarized by its 
FWHM, as it depends on the structure function of the phase 
residuals, and in particular on the compensation order. The 
more important this latter, the more intense is the "coher- 
ent core" (with a width of X/D) on top of the PSF halo. 
Deconvolution processes can then be applied, whose perfor- 
mance increase together with signal-to-noise, i.e. with the 
restoration quality. Expecting this high frequency restora- 
tion makes it mandatory to obey the Nyquist critera for the 
image sampling. 

Expectations are different when studying faint galaxies 
with 3D spectroscopy at i? « 10000. The intrinsic faintness 
of the dispersed signal prohibits, anyway, any attempt to 
spatially sample it at Nyquist frequency. Hopefully for us, 
this latter is not an issue, since we are not interested in 
the ultimate telescope resolution, but only in a moderate 
resolution of 0.25 arcsec. Our concern is therefore to max- 
imise the flux from the object within a square aperture of 
0.25 X 0.25 arcsec? -our spatial resolution element. This has 
two consequences. Firstly, it increases the signal from the ob- 
ject within the spatial resolution element, whereas the noise 
level (especially thermal and sky background) will remain 
constant: the spectroscopic SNR will increase. Secondly, the 
flux spread in the neighbouring pixels will be minimised, 
thus avoiding spatial contamination. 

These improvements of the image quality expected in 
spectroscopy (resolution, SNR improvement, reduction of 
confusion) are diflicult to relate with the parameters com- 
monly used in adaptive optics to characterise the quality of 
the restored image. Classical parameters are, as an exam- 
ple, the full width at half maximum (FWHM) or the Strehl 
ratio (SR ; this latter is deflned as the central intensity of a 
restored PSF, divided by the central one of the diffraction- 
limited PSF, taking both PSFs with the same total energy). 
For 3D spectroscopy applications, the criterium tradition- 
ally used, and that defines the spectroscopic image quality, 
is in fact the fraction of the total PSF flux ensquared within 
the spatial resolution element. We will therefore use this 
criterium, the ensquared energy, in the following sections of 
this paper. 



Zernike mode index 
2 16 36 



case A 



case B 



Figure 1. The hatched area of the bars indicate the spatial 
frequency range of the residuals, described by the index of the 
Zernike mode. Case A has been chosen because it mimics a case 
typical of a low-order AO system. Numbers of case B (residuals 
is the range Ziy—sg) are chosen to provide the same wavefront 
flatness (same phase variance) as in case A. The two PSF will 
have the same phase variance, same SR, and same FWHM. 



3.2 The need for a high-order compensation. 

FALCON had initially began with the wrong idea that a 
low-order, or a moderate compensation should be sufficient 
to achieve a goal that, at least at the first glance, appeared 
to be modest : obtain a resolution of 0.25 arcsec. We demon- 
strate in the following sections that getting such a resolution 
implies a finer spatial sampling (0.125 arcsec) compared to 
existing instruments, and consequently a huge loss in signal- 
to-noise per spatial sample. This loss has to be compensated 
by an increase of the ensquared energy. We point out that a 
high ensquared energy can be obtain only by using a high- 
order compensation. This effect is implicitly demonstrated 
at the end, in the last result section of our article, where 
we will see that the ensquared energy still take advantage 
of a high-order correction. However, we want to insist on 
this particular point from now, because this explains why 
we have especially focused a large part of our study on the 
number of freedom of the system shall have, and why we 
study all the compensation levels up to 120 Zernike modes 
in this paper. 

We illustrate here how the ensquared energy is a pa- 
rameter sensitive to the spatial frequency range of the phase 
residuals. We will compare two illustrative cases : they are 
comparable, in the sense that their wavefront fiatness is the 
same (i.e. same amount of residual phase variance). The two 
cases just differ by the spatial frequency of the phase resid- 
uals. 

Any publication about AO compensation focuses on 
phase residual variance, which is considered as the key pa- 
rameter that optimizes the image quality. We exhibit here 
two striking examples of PSFs where the ensquared energy 
is fundamentally different (more than a factor of 2), while 
their associated phase variance, Strehl ratio, and FWHM 
are the same. 

For our illustration, we performed some Monte-Carlo 
numerical simulations of an AO corrected PSF, where the 
wavefront (from a fully developped Kolmogorov spectrum) 
is corrected by zeroin g some coef ficients of its expansion into 
Zernike polynomials l|Nol]|[l976h . As we just study the in- 
fiuence of the compensation order, here we focus only on 
the spatial aspect of the correction, ignoring errors due 
to time delay, anisoplanatism or noise. The conditions of 
the simulation are a 8-meter telescope (VLT case), a see- 
ing of 0.81 a r csec a t 0.5 /im (median seeing of VLT site, 
iMartin et al.l l|2000l )). leading to a ro equal to 52 cm at 
the imaging wavelength of 1.65 /im (H band). We com- 
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AO-corrected PSF (low order correction) 



• 



AO-corrected PSF (high order 




Figure 2. AO-corrected PSFs in the case of low order correction 
(perfect correction of ZPs up to jmax = 36, left) and high order 
correction (perfect correction of ZPs 2 to 16 and 37 to oo, right). 
8 meter diameter telescope, seeing of 0.81 arcsec at A = 0.5 fim, 
H band (A = 1.65 fim). The simulated field of view is equal to 
1.70 X 1.70 arcsec^. 



PSF radial profiles 



Ensquared energy protiie; 




Distance from center (arcsec) 



Aperture size (arcseo) 



Figure 3. Radial profiles (left) and ensquared energy profiles 
(right) of the AO-corrected PSFs. The solid line corresponds to 
low order correction, the dashed line to high order correction. 



pared two clear-cut cases of correction, similar since they 
both lead to the same amount of residual phase variance 
ffres = 0.0132(D/ro)^'''' rad^, but different because the com- 
pensation concerns different domains of spatial frequencies, 
as illustrated in figure [T] 

• for the case A, we mimic a typical low-order AO system, 
and we compensate for the polynomials Z2 to Z^e. Phase 
residuals are high-order perturbations Z>37. 

• for the case B, we consider the phase residuals are 
now low-order perturbations in the range Zn-^Q. Those last 
numbers are chosen so that the spatial frequency contents 
of the phase does not overlap with case A, but has exactly 
the same variance. 

The AO corrected PSFs in H band are shown on the 
figure [21 whereas figure [3] shows the radial (left) and en- 
squared energy profiles (right) . Those figures aim to demon- 
strate how different the shapes of the PSFs can be, al- 
though they simultaneously exhibit identical residual phase 
variances, identical central intensities (predictable from the 
phase variance: both Strehl ratios are close to 30%), and 
identical FWHMs: the shapes of the coherent cores are iden- 
tical, only the energy distribution of the halo is different, 
because it is strongly determined by the spatial contents of 
the phase residuals. 

In terms of astrophysical implications, a PSF such as in 
case B (low-order residuals, right of figure[2]) would be catas- 
trophic for exoplanet detection, where the intensity in the 
first rings has to be minimised because of contrast issues. On 
the opposite, such a PSF is interesting for 3D spectroscopy 



as the ensquared energy is increased: this is particularly no- 
ticeable on the right of the figure |3] 

One has therefore to keep in mind that for a fixed 
amount of phase variance removed, the ensquared energy 
will strongly be enhanced in the case of a high order 
compensation, which is particularly effective at reducing 
the halo compared to low-order modes correction. The 
conclusion is that the common rule of thumb, according 
to which the phase variance error budget should be bal- 
anced between undermodeling error and other sources, 
becomes obsolete when it is a matter of designing an AO 
system for spectroscopy. Despite the moderate resolution 
(« 0.25 arcsec) required, a rather large number of actuators 
is required to bring the fight from the further parts of the 
halo into the inner parts, so as to increase the spectroscopic 
SNR, and only a high-order modes compensation can allow 
us to achieve this. 

However, compensating for high order modes is not go- 
ing without some difficulties from the AO point of view, the 
major limitation being anisoplanatism: high order modes 
unf ortunately decorrelate faster with the separation an- 
gle (|Chassatlll989l : lRoddier et al.lll993l ). Unfortunately, the 
usual observing conditions of distant galaxies are at high 
galactic latitudes, where natural guide stars become scarce 
and are unlikely to lie within the isoplanatic patch ^0. This 
requires to work most of the time with large separation an- 
gles, that precisely forbid the compensation of high-order 
modes. This drawback makes classical AO unappropriate 
for the statistical studies of distant galaxies. 

In the section [J] we will describe a new approach of 
adaptive optics, which solves anisoplanatism problems, 
which allows to use natural guide stars at distances greater 
than 80, and which is still efficient to correct high order 
modes. 

Before this, we want to be able to establish requirements 
about the ensquared energy needed for a proper statistical 
study on high redshift galaxies, in order to derive later some 
specifications for the AO system we propose. 



3.3 Case study: 3D spectroscopy of UGC 6778 
observed at high redshift 

Thanks to numerical simulations, we show in the following 
subsections some quantitative results about the gain in 
spectroscopic SNR that can be achieved when AO and 3D 
spectroscopy are used together. We would like to emphasize 
here that we do not attempt to perform a complete simula- 
tion of a 3D spectrograph feeded by an AO system, because 
of the lack of 3D spectroscopic data currently available in 
the litterature for high redshift galaxies, and especially in 
the near-infrared wavelength range, where the Ha emission 
line is observed for redshifts z ^ 0.5. We therefore do very 
simple assumptions about the instrument as well as the 
observed objects, in order to compute first estimates about 
the SNR achieved thanks to AO correction. Such results 
can be used to derive very preliminary required ensquared 
energy values, that can then be translated into instrumental 
requirements for the AO system. A ACDM cosmology with 
Ho — 70 km/s/Mpc, i},n = 0.3 and JIa = 0.7 is assumed in 
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=1 from lGarrido et al.l ll2002h 
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from ISaggett et al. I l ll998t) 



Table 1. Photometric and dynamic characteristics of UGC 6778. From left to right: the identifier of the 
galaxy, its morphological type, the recession velocity corrected from the local's group infall toward Virgo 
VviR, the distance D in Mpc, the absolute magnitude in B band Mb, the exponential disk scale lenght 
i? D , the disk surface brightness /x/, , the inclination i and the position angle PA, and the maximum velocity 
observed on the edge of the disk 



the rest of this paper. 

We used for this study the Ha image of the UGC 
6778 galaxy provide d by member s of th e GHASP survey 
(|Garrido et al.l I2OO2I , |2003| . |2004 |2005| ). with a spatial 
sampling so = 0.960 arcsec/pixel. The table 1 lists the 
photometric and dynamic properties of this galaxy. These 
data were then used to simulate two science cases, where 
we assumed that the galaxy would be respectively located 
at z = 0.9 OT z = 1.5, and that a microlens+fibers integral 
field spectrograph similar to GIRAFFE would be used to 
perform 3D spectroscopy on this galaxy. 

The first steps required to perform these simulations are 
the following: 

(i) Use of the original Ha image to simulate the distant 
galaxy image, taking into account the decrease of the appar- 
ent diameter due to redshift as well as the change of spatial 
sampling (see iGiavalisco et al. I ll996l for more information). 
We assumed a spatial sampling in the simulated galaxy im- 
age equal to half of the telescope's diffraction limit (Nyquist 
sampling), i.e. X/2D. 

(ii) Computation of the apparent magnitude m in the ob- 
serving band (J or H) from nearby galaxy's absolute mag- 
nitude M, using distance modulus, and taking into account 
the morphological K-correction K{z) due to redshift into the 
distance modulus equation: m — M — 51og(DL/10) +K{z), 
where Dl is the lumi nosity distance. We used values from 
iMannucci et al.l l|200ll ) for the K-correction. The flux in the 
continuum fcont (in ergs s~^cm~^A~^) can then be com- 
puted from the apparent magnitude m. 

(iii) Computation of the total flux (in ergs s^^cm"^) of 
the galaxy into the redshifted Ha line, and normalisation of 
the simulated Ha image. 

(iv) Simulation of the distant galaxy continuum image 
assuming pure Freeman exponential profile, taking into ac- 
count the disk scale length Rd and the angular distance Da- 
The flux normalisation is done using fcont- 

The figure |4] shows therefore the original Ha image 
coming from the GHASP survey (z ~ 0.0, pixel size 
So = 0.96 arcsec), as well as the simulated Ha and contin- 
uum map corresponding to the first science case {z = 0.9, 
pixel size s^ = 0.016 arcsec) - 

We are going now to focus on the two HI I regio ns on 
the bottom left of the Ha image. ICarrido et al.l (|2002l ) have 





Original Ha image 
z RJ 0.0, A = 0.656 /im 
2.62 X 2.62 arcmiri^ 
So = 0.96 arcsec 



Simulated Ha image Continuum map 
2 = 0.9, A = 1.250 Atm 
1.80 X 1.80 arcsec^ 
Sz = 0.016 arcsec 



Figure 4. Original Ha image (left), simulated Ha (middle) 
and continuum map (right) of the UGC 6778 galaxy as if it was 
observed at 2 = 0.9 (in J band), with the same physical size than 
today. The simulated field is 1.80 X 1.80 arcsec? - 



indeed shown that they have very similar velocities (~ 890 
km/s). We therefore do the assumption (sufficient for the 
rest of this study) that the integral field spectrograph we 
use has a sufficient spectral resolution K so that the flux 
sampled around the wavelength Ao ( 1 -I- w/c) (where Ao is the 
observing wavelength, i.e. 1.25 or 1.65 /im) is the flux only 
emitted in the Ha line by these two regions, v being their 
average velocity. We therefore isolated these two regions in 
the simulated images we created. 

However the simulated Ha and continuum images sup- 
pose infinite spatial resolution, and do not take into account 
atmospheric seeing or AO partial correction effects. 

We therefore convolved those images with different 
point spread functions (PSF) corresponding to all the pos- 
sible orders of AO correction, from to 120 Zernike, at 
1.25 \im and 1.65 /.im. The principle of the simulation is 
the same than in section [3^ and we assumed the same con- 
ditions (8m telescope, seeing of 0.81 arcsec at A = 0.5/^m, no 
anisoplanatism, noise or temporal error). As a consequence 
of this AO compensation, a better contrast in the images 
as well as a better ensquared energy is observed, as already 
explained in section [3Tl1 The result of the AO impact is il- 
lustrated, for a few cases, on figure [5] 

Those images could then be used to compute the flux 
entering into the microlens sampling each HII region, and 
then compute the final number of electrons on each pixel of 
the detector located at the output of the spectrograph, re- 
quiring therefore to know the spectral proflle of the line. We 
assumed a gaussian line profile for the line, being the con- 
sequence of some internal velocity dispersion ay in those 
HII regions, as well as the finite spectral resolution K of 
the spectrograph. As a result, the line profile on the spec- 
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Figure 5. Illustration of the effect of adaptive optics com- 
pensation on the images. The simulated images cover a field of 
0.50 X 0.50 arcsec^ (4x4 microlenses of 0.125 X 0.125 arcsec^). 
A analysis has been done, that covers compensation orders from 
to 120 Zernike modes. A few cases have been picked up and 
are presented here for illustration purpose. From left to right; 
diffraction-limited image, seeing-limitod image, AO-corrected im- 
age for 46 and 100 ZPs. 



trograph detector has a gaussian shape, with a standard 
deviation atot = \J'^v,\ + ""^.^ — Xo x o-y/c being the 
standard dispersion in wavelength due to the internal veloc- 
ity dispersion of the HII regions, and as = Ao/2.354i? being 
the standard deviation due to the spectral resolution, with 
Ao being the observing wavelength and c the light celerity. 
Therefore, if we call Fml the total flux in HII region sam- 
pled by the spectrograph's microlens, the line's profile /(A) 
as a function of wavelength will be written: 



fW = 7^ exp - 



(A-Ao)^ 



(1) 

This line profile is then going to spread over two direc- 
tions on the spectrograph's detector: the horizontal direction 
will correspond to the spectral information, and the perpen- 
dicular direction will correspond to the spatial information. 
This allows us to compute the number of electrons due to 
the Ha line per pixel on the spectrograph detector, and 
then the spectroscopic SNR. If AA is the spectral sampling 
of the detector, Stot the telescope area, iexp the exposure 
time, and Ttot the total transmission of the system (which 
takes into account the detector efficiency), the number of 
electrons Nna coming from the redshifted Ha line on the 
pixel sampling the line at the Ao wavelength is then equal 
to: 

Nhc, = 2 / ./■(A)rfA X Stot X t^xp X Ttot x Xo/hc (2) 
Jo 

where h and c are respectively the Planck constant and 
the speed of light. 

We then have the fiux from the line going onto the de- 
tector pixel. However, as we are interested in computing 
the spectroscopic SNR, we therefore need to take into ac- 
count all the sources of noise, which can be divided in three 
categories: the line and continuum photon noise, the back- 
ground photon noise, and the detector noise. As both line 
and continuum signals have poissonnian statistics, therefore 
their associated photon noise has a variance equal to their 
respective value, i.e. Ns and N^ont, where the number of 
electrons Ncont due to the continuum is equal to: 



iVc, 



fcont X Stot X te.:cp X AA X Ttot X A0//1C (3) 



fcont being the flux from the continuum entering the 
microlens after convolution of the continuum image by 
the AO corrected PSF. Concerning the background noise, 
this latest is the sum of the sky background as well as the 
instrumental thermal background. Then, if we call Nb the 




Figure 6. Explanation of the effective flux notion. The first 
image on the left shows the original image, with the two HII 
regions and the superposed microlens (in dashed line) . The second 
image shows only the area sampled by the microlens. The third 
and fourth show this latter convolved by the uncorrected PSF 
(third) and by an AO-corrected PSF (fourth image, correction of 
60 ZPs). The effective fiux from the Ha line is then equal to the 
fiux of the region after convolution by the PSF going into the 
microlens. 



total number of electrons coming from the backgroud, as 
it has also Poisson statistics, its associated photon noise 
will then have a variance equal to Nb- At last, detector 
noises, i.e. dark current noise and readout noise, have to be 
considered. 

Knowing all these values, we can then compute the SNR 
per spatial resolution element and per spectral pixel, which 
will be written: 



SNR- 



^/Nhc, + Noont +Nb+ Np,x.spat{ND + iV.:ipCr|) 

(4) 

with: 

• Nb the number of electrons coming from the instru- 
mental and sky background 

• Nd ~ DC X Texp the number of electrons dues to the 
dark current DC (in e~ /s/pix) 

• Nexp the number of elementary exposures 

• (T_H the readout noise (in e~ /pix) 

• N pix, spat the number of pixels on which the fiux is 
spread in the spatial direction on the detector 

The equation |4] therefore gives the expression SNR on 
the measured quantities. However, the figure [5] shows that 
as the AO correction improves, the image contrast becomes 
better, with an increase of the effective light coming from 
the HII region. As a result, the flux going into the microlens 
sampling each HII region is therefore the sum of two terms: 
the effective flux emitted by the HII region, and a pollu- 
tion term coming from adjacent regions, which is going to 
degrade the spectrum of the area sampled by the microlens, 
provided as an output of the 3D spectrograph. This has 
leaded us to define another signal-to-noise ratio, which we 
will call Effecttve signal-to-notse ratio (ESNR), and which 
gives a measure of the effective gain provided by the coupling 
of AO with 3D spectroscopy. Its expression is: 

ESNR = jVg^.e// 

^TVhc + A^cont +Nf+ Np,^,spat{Nn + N^^pol) 

(5) 

where NHa,eff is the number of electrons due to the 
effective Ha fiux emitted by the area sampled by the 
microlens after convolution by the AO corrected PSF, 
with no contribution from the adjacent areas around the 
microlens. This notion of effective flux is explained on the 
figure m 
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One will notice that denominators in the equations |4] 
and [5] are the same, and therefore that the ESNR is a lower 
bound of the SNR. This is explained by the fact that in the 
case of the ESNR, all the noise sources must appear in the 
denominator. Concerning the line, this latter is perturbed 
by a noise, which is the sum of two terms: the photon noise 
of the effective flux, and a pollution term due to all the light 
coming from the area around the microlens and entering the 
microlens because of the convolution. As both quantities 
have poissonnian statistics, and as the convolution operator 
is distributive, the sum of these two terms is equal to the 
photon noise of the measured light entering the microlens. 

We would however like to emphasize here that the 
ESNR defined at the equation [5] and the classical SNR 
definition shown in equation |4] cannot be directly compared, 
in particular when the ensquared energy changes. As it 
is shown in figure [SJ the ESNR numerator is the flux 
coming only from the region sampled by the microlens 
after convolution by the AO corrected PSF, whereas the 
SNR numerator is the flux not only coming the microlens, 
but also the flux from adjacent areas and polluting the 
microlens because of the convolution. As different images 
are convolved for ESNR and SNR computations, there is no 
direct relation (i.e. multiplication) between these numbers. 
Moreover, as explained in the section 13.11 an increase of 
the ensquared energy is due to a better compensation order 
provided by the AO system. Thus different ensquared en- 
ergy values correspond to different AO corrected PSFs, and 
the results of the convolution of the same high-resolution 
galaxy image by these different PSFs also cannot directly 
be comparable. 

Once these two different quantities (SNR and ESNR) 
defined, we performed some numerical simulations in order 
to quantify the gain in ESNR provided by the combination 
of AO with 3D spectroscopy. The results of these simulations 
are shown in the next section. 

3.4 Results 

Thanks to simulations, we are now able to show some 
preliminary results about the improvement in ESNR thanks 
to AO. Two science cases, corresponding to the observation 
of UGC 6778 at two different redshifts, are considered in 
the two next subsections. 

For the first one, we assumed that the observed galaxy 
would be located at z = 0.9. In that case, the Ha emission 
line is redshifted at a wavelength of 1.25 ^m, corresponding 
to the central wavelength of J band. Moreover we assumed 
that the galaxy had the same properties than the distant 
large spirals obs erved in the CFRS l|Lillv et al.1 1 19981 : 
IZheng et al.l |2004| ) . In particular the physical size of the 
simulated galaxy is the same than the one observed in 
the local universe. As there are no direct measurements of 
Ha fiuxes available in the litterature for CFRS galaxies at 
these redshift, we there fore computed the H a luminosity 
from [O/J] 3727 fiuxes l|Hammer et al.1 1 19971 '). in order to 
normalise the simulated Ha image. 

For the second one, we assumed that the observed 



galaxy would be located ai z — 1.5. In that case, the Ha 
emission line is redshifted at a wavelength of 1.65 /im, 
corresponding to the central wavelength of H band. More- 
over we assumed that the galaxy had Ha fiuxes similar 
to the ones ob served on the popula tion of " 'BM galaxies" ' 
introduced by ISteidel et al.l l|2004 ). The apparent size of 
the galaxy was then chosen to cover a field of 1 arcsec^, 
consistent with apparent sizes obs erved on distant gala xies 
at these redshifts in deep surveys (|Bouwens et al.|[2004l ). 

We focus now on the two HII regions we used to 
compute the SNR and ESNR. Together they contribute at 
10% of the total Ha flux in the simulated Ha image. We 
also assumed that they had an internal velocity dispersion 
with FWHMv = 30 km.s"^ (eg cry ^ 13 km.s"^), leading 
as explained previously to a dispersion in wavelength on 
the spectrograph's detector equal to (ta,v = Aq x cry/c, Ao 
being equal to 1.25 or 1.65 /im. The spectral resolution was 
then flxed to _R = 10000, allowing to resolve the velocity 
dispersions of the HII regions. This then allowed us to 
deflne two cases of spectral sampling, corresponding to 
the horizontal direction on the spectrograph's detector: 
AA = 0.3 A/pix (for J band) or AA = 0.6 A/pix (for 
H band). Such spectral samplings allow us to sample 
the spectral resolution with more than two pixels, thus 
providing a better sampling than Nyquist sampling, and 
then do a precise estimation of the barycenter of the Ha 
emission line in order to recover the galaxy's velocity field. 
For the vertical direction on the spectrograph's detector, 
corresponding to the spatial information, we considered a 
number of pixels Npix,spat ~ 3, similar to the number of 
pixels used on the GIRAFFE spectrograph's detector. Such 
a number of pixels, here again better than Nyquist sam- 
pling, is quite convenient to resolve spectrum deblending 
problems on the spectrograph's detector. 

For the sky background, we assumed that the spec- 
trograph is observing between atmospheric OH lines, 
so the sky noise is equal to the sky continuum. ICubvl 
(2000) has measured the sky continuum at the Cerro 
Paranal observatory, and found respective values of 1200 
ph/ s/m^ / fj.m/arcsec'^ at a wavelength of 1.25 fj,m, and 
2300 ph/ s/m^ / fj.m/arcsec'^ at a wavelength of 1.65 fj,m. We 
used the characteristics of the VLT's ISAAC spectrograph 
to simulate the instrument: the instrumental background 
was assimilated to a blackbody of temperature T — 273K 
and an emissivity equal to 25%, a total transmission 
(atmosphere-|-telescope-|-instrument-|-detector) of 13% at 
1.25 jj,m (central wavelength of J band) and 19% at 1.65/im 
(central wavelength of H band) was assumed, and the 
detector had a dark current DC = 30e — /h/pix and a 
readout noise RON = 10 e — /pix. 



3.4-1 First science case: 3D spectroscopy of UGC 6778 at 
z^O.9 

Here we give the results of simulations of the 3D spec- 
troscopy of UGC 6778 as if it was observed at 2; = 0.9. 
We assumed that the galaxy had the same physical size as 
it has now. However the fiux in the Ha line is consistent 
with the one observed of distant galaxies. Such physical 
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conditions correspond to the ones observed on large spiral s 
in the distant universe (|Lillv et al.lll998l : IZheng et al.ll2004 ). 

The original Ha image covers a field of 2.62 x 
2.62 arcmin^, and the galaxy has a distance D — 17.04 
Mpc, which leads to a physical field of di x di kpc^, with 
di = 13.40 kpc. Now, if we assume that the same galaxy has 
a redshift z, we can then compute its luminosity-distance Dl 
and its angular distance Da, and we found that the observed 
field has therefore an apparent size equal to 6'z x 9^ arcsec? , 
with 9z = d\/DA = 1.80 arcsec. Moreover, the spatial 
sampling is not the same as we simulate the galaxy as if 
it was observed by the VLT (8 meter diameter telescope) 
at its diffraction limit, meaning the pixel sampling Sz is 
equal to Nyquist sampling (A/2_D), leading to a pixel size 
Sz = 0.016 arcsec. Therefore the number of pixels in the 
original and simulated images are not the same: as shown 
in figure 131 the original image has 164 x 164 pixels, whereas 
the simulated image has 111 x 111 pixels. 

This simulated image must be then normalized so that 
its flux is equal to the flux in the Ha line. As explained 
before, there are not direct measurements of Ha fluxes on 
the CFRS galaxies, theref ore we used [O J J] 372 7 luminosi- 
ties to compute Ha fluxes. iHammer et al. I l|l997l ) found that 
the average rest-frame equivalent width of the [077J3727 
line in CFRS galaxies located at z = 0.9 was equal to 
£;W([0//]3727) = 40A. Moreover. iKennicuttI (|1992 ) found 
that there is a relation between the equivalent width in the 
[07/]3727 and Ha line, in the form EW{Ha) = 2.5 x 
EW{[OII\'i727). We therefore assumed that the Ha line 
had a rest frame equivalent width EW{Ha) = lOOA. 

Once EW{Ha) is known, the flux in the contin- 
uum fcont must be computed to know the flux fna 
in the Ha line, as both are linked to the relation: 
fHa = .fcont X EW {H a) X {I + z) . fcont IS dlrcctly linked 
to the apparent magnitude of the galaxy. We see on the 
table 13.31 that UGC 6778 has an absolute magnitude 
in B band Mb = —20.6. However, as the Ha line is 
redshifted in J band, we are interested by the galaxy's 
absolu te magnitude in J band. For a SBc galaxy like UGC 
6778, |m annucci et all liool) give the following colours: 
{B - V) = 0.70, {V - K) = 3.03, {J - H) = 0.66 et 
{H — K) = 0.25. We t hen find an abso l ute re st-frame mag- 
nitude M,j = -23.42. iMannucci et all (120011 ) give also the 
following K-correction for a Sc galaxy: K{z = 0.9) = 0.155. 
The redshift z can then be used to compute the luminosity 
distance Dl and the distance modulus. We therefore find 
an apparent magnitude J = 20.41, corresponding to a flux 
in the continuum fcont ~ 8.09 x W~^^ ergs/crn^ / s/ A, 
leading to a Ha flux of 4.17 x 10~^^ erg/s/cm?. This flux 
was then used to normalise the simulated Ha image. 

As explained before, the continuum emission has also 
to be taken into account. This latter is not equally spread 
over the galaxy, as the intensity of dis k in spiral gala xies 
follows a decreasing exponential law l|Freemanl Il970l ). in 
the form 7(r) oc exp"'"''^'', Rd being the disk scale length, 
whose value for UGC 6778 is given in table 13.31 It is 
therefore possible to compute Rd,z, the angular size of 
of the disk scale length when the galaxy is observed at 
a redshift z, by the relation Rd,z = Rd/Da, Da being 
the angular diameter distance, and we found a value of 
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Figure 7. Measured spectroscopic SNR (top) and effective spec- 
troscopic SNR (ESNR, bottom) on the Ha line for a microlens 
of 0.125 arcsec as a function of the ensquared energy in a 
0.25 X 0.25 arcsec^ square aperture. The considered galaxy is 
located at z = 0.9 (observation in J band) and has the same char- 
acteristic of distant large spirals. Each curve corresponds to the 
average of the SNR of each microlens sampling the HII regions, 
microlens, with exposure times going from 1 hour to 8 hours. The 
dashed line corresponds to SNR = 3. 



Rd,z = 0.25 arcsec. This value, as well as the position 
angle PA of the galaxy and the inclination i (which are 
given in the table 1), were used to compute the continuum 
map of the galaxy. This latter was then normalised so that 
its total flux was equal to fcont- The simulated Ha image 
as well as the continuum map are shown on the flgure]?] 

Once those two images created, it is possible to compute 
the measured spectroscopic and effective SNR per spectral 
pixel. However, as our goal is to show the gain given by 
the combination of AO with 3D spectroscopy, we convolved 
these latter images by simulated AO corrected PSFs, with 
a degree of correction going from ZPs (seeing-limited 
PSF) to 120 ZPs, leading to an increase of the ensquared 
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energy as the number of corrected modes increases. We 
also considered exposure times between 1 and 8 hours, with 
individual exposure times of 1 hour. The figure [7] shows 
the evolution of these SNRs as a function of the ensquared 
energy into a 0.25 x 0.25 arcsec? square aperture. More 
precisely, if we look at the evolution of the ESNR which tn 
fine says the real improvment provided by AO correction, 
we see that a minimum ensquared energy of 26% allows 
to reach an ESNR of 3 after an exposure time of 8 hours. 
However, a slightly increase to an ensquared energy of 30% 
allows to reach the same ESNR value, but only after 3 
hours of exposure time. 

As it was shown previously, the AO correction not only 
increases the ensquared energy into the spatial resolution 
element, but it also increases the sharpness of the AO 
corrected PSF, therefore providing a better image contrast. 
This is indeed shown on the figure |8l where it can be seen 
that the two HII regions start to be well separated after 
the perfect correction of the first 46 Zernike Polynomials. 
This level of AO correction corresponds to an ensquared en- 
ergy equal to 28% into a 0.25 x 0.25 arcse(? square aperture. 

This study has therefore shown that the combination 
of 3D spectroscopy with AO is definitely effective in terms 
of angular resolution and SNR improvement. This means 
that the 3D spectroscopy of 2 ~ 1 galaxies with an angular 
resolution of 0.25 arcsec (definitely better than atmospheric 
seeing) should be achievable on a 8 meter telescope located 
on the ground such as the ESO VLT. This requires the AO 
correction to provide a minimum ensquared energy of 28% 
in a 0.25 x 0.25 arcsec^ square aperture. In that case, a 
spectral resolution R = 10000, allowing to resolve velocity 
dispersions with FWHMv ^ 30 km/s, can be used, as an 
effective spectroscopic SNR of 3 can be reached after an 
exposure time of only three hours. 



3.4-2 Second science case: 3D spectroscopy of UGC 6778 
at z — 1.5 

We give here the results of simulations of the 3D spec- 
troscopy of UGC 6778 as if it was observed sd, z — 1.5. We 
assumed here that the galaxy covers an apparent surface of 
1 arcsec^, which is consistent with apparent sizes measured 
on deep surveys. This apparent surface was then used to sim- 
ulate the apparent Ha and continuum images. We are then 
going to focus on the continuum and Ha flux normalisation. 

We used the data from lErb et all l|2003h and 
ISteidel et al.l l|2004l ). who observed with the LRIS spectro- 
graph on the Keck telescope a sample o f BM/BX galax - 
ies located between 2 = 1.4 and z = 2.5. lErb et al.] (120031 ) 
found on their sample {{z) — 2.28) an average apparent 
magnitude (R) — 24.37, and an average apparent Ha flux 
f{Ha) = 4.6 X 10"^'' ergs cm~^. One wiU notice that 
the apparent Ha flux is comparable to the one computed in 
the previous science case [z — 0.9) where we simulated the 
observation of a distant large spiral galaxy belonging to the 
CFRS. In fact the two science cases have been chosen using 
previous well-known studies of distant galaxies. At 2; = 0.9, 
the CFRS has studied average L* galaxies, while at z = 1.5, 



lErb et all l|2003l ) have focused on a sample of specifically 

Ha luminous g alaxies. 

Moreover ISteidel et al] l|2004l ) found that galaxies 
with a redshift (z) = 2.23 ± 0.31 had an average color 
{R — K) = 3.25 ±0.53. We therefore found that those galax- 
ies had an apparent magnitude {K = 21.12). It is therefore 
possible to compute their K absolute magnitude Mk 
using again the distance modulus, assuming a k-correction 
K{z) — —0.57 at 2 = 2.3 (extrap olation to 2 = 2 .3 of the 
k-correction in K band given by iMannucci et al.l for a Sc 
galaxy), leading to Mk ~ —24.3 8. Moreover, as Sc gal axies 
have a color {H - K) = 0.25 (|Mannucci et all 120011 ). we 
found that such a galaxy has a rest-frame H absolute 
magnitude equal to Mh ~ —24.13. Assuming a k-correction 
K(z) = 0.078 at 2 = 1. 5 for a Sc galaxy observed in H band 
(^Mannuc ci et al.l l200lf ). we therefore found an apparent 
magnitude H = 20.99 for the case where such galaxies 
would be located ed z — 1.5, corresponding to a fiux in the 
continuum fcont ~ 1-67 x 10~^^ erg/s/cm^ /A. 

These flux values were used to normalise the simulated 
images of the distant galaxy. Then we did the same 
computations as previously explained, in order to compute 
the measured spectroscopic and effective SNR per spectral 
pixel on the spectrograph's detector, and especially show 
the improvements provided by AO correction. This time we 
considered exposure times from 3 to 24 hours. The result of 
these simulations is shown on figure 1101 where we can see 
that an effective SNR of 3 can be reached after a minimum 
exposure time of 9 hours, requiring an ensquared energy 
after AO correction equal to 60% into a 0.25 x 0.25 arcsec'^ 
square aperture. The same ESNR value can be reached 
with an ensquared energy after AO correction equal to 
40%, requiring then an exposure time of 18 hours, or also 
with an ensquared energy after AO correction equal to 
33%, requiring then an exposure time of 24 hours. However, 
in the previous section we saw that a minimum ensquared 
energy must also be reached so that the two HII regions can 
be separated. This is explained in the figure ^ where we 
can see that the two HII regions start to be separated for 
an ensquared energy after AO correction better than 30%. 

This study has therefore shown that the combination 
of 3D spectroscopy with AO should allow also to perform 
the 3D spectroscopy of galaxies located at 2 = 1.5, also 
with an angular resolution of 0.25 arcsec (definitely better 
than atmospheric seeing) and with a spectral resolution R = 
10000, allowing to resolve velocity dispersions FWHMv ^ 
30 km/s, as well as reaching effective spectroscopic SNR of 
3. This can be achieved with a minimum ensquared energy 
after AO correction of 35% and an exposure time of 24 hours, 
knowing that the same ESNR value can be reached with an 
ensquared energy of 40% and a shorter exposure time (18 
hours). 

3.5 Conclusion 

We studied in this section the improvement provided by 
the combination of 3D spectroscopy with Adaptive Optics 
(AO). Firstly we have shown that thanks to AO, it is 
possible to reach angular resolution better than seeing, 
but also to increase the spectroscopic signal-to-noise ratio 





Figure 8. Images of the HII regions convolved witli tlic AO corrected PSF, for different levels of AO correction in J band. The title of 
each image sliows tlic number of corrected modes and tiic corresponding ensqucired energy into a 0.25 X 0.25 arcse(? . The field of view 
is equal to 0.50 x 0.50 arcse(? (4 x 4 0.125 arcsec microlenses) . 
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Diff. limited image 



Jmax=0, ensquared energy (0.25 arcsec) = 1 7.4 % Jmax=5, ensquared energy (0.25 arcsec) = 21 .2 % 




Jmax^l 3, ensquared energy (0.25 arcsec) ^ 26.9 % Jmax=1 9, ensquared energy (0.25 arcsec) = 29.7 % Jmax=29, ensquared energy (0.25 arcsec) ^32.5 % 




Jmax=37, ensquared energy (0.25 arcsec) = 34.9 % Jmax^46, ensquared energy (0.25 arcsec) = 38.0 % Jmax^67. ensquared energy (0.25 arcsec) ^ 46.2 % 




Jmax=85, ensquared energy (0.25 arcsec) = 51 .9 % Jmax=99, ensquared energy (0.25 arcsec) = 55.5 % jmax=1 20, ensquared energy (0.25 arcsec) = 60.2 % 



Figure 9. Images of the HII regions convolved witii tlie AO corrected PSF, for different levels of AO correction in H band. The title of 
each image shows the number of corrected modes and the corresponding ensquared energy into a 0.25 X 0.25 arcsec^ square aperture. 
The field of view is equal to 0.50 x 0.50 arcsec^ (4 x 4 0.125 arcsec microlenses) . 
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Figure 10. Measured spectroscopic SNR (top) and eflfective 
spectroscopic SNR (ESNR, bottom) on the Ha line for a mi- 
crolens of 0.125 arcsec as a function of the ensquared energy in a 
0.25 X 0.25 arcsec^ square aperture. The considered galaxy is lo- 
cated at z = 1.5 (observation in H band) and covers an apparent 
field of 1 arcsec^ . Each curve corresponds to the average of the 
SNR of each microlens sampling the Hll regions, with exposure 
times going from 3 hour to 24 hours. The dashed line corresponds 
to SNR = 3. 



of 3 can be reached for a galaxy located a,t z = 0.9, requiring 
a minimum ensquared energy into a 0.25 x 0.25 arcsec^ 
after AO correction equal to 30% and an exposure time of 
3 hours. The same ESNR value can be reached for a galaxy 
located at z — 1.5, requiring then a minimum ensquared 
energy after AO correction of 35% and an exposure time of 
24 hours, however a slight increase of the ensquared energy 
to 40% should allow to reach the same ESNR value after 
an exposure time equal to 18 hours. 

A study of the AO corrected PSF has shown that high 
order modes need to be corrected to improve the ensquared 
energy. This has been confirmed by our simulations of the 
3D spectroscopy of distant galaxies, where we found that 
at least 46 Zernike polynomials need to be corrected to 
reach an ensquared energy of 30% in J band and 40% in 
H band (such results are not definitive, as we used the 
dynamical properties of one galaxy observed in the local 
universe as well as those of a few distant galaxies, but they 
can be used as a first basis, the data required to perform 
such studies being still rare in the community). Such a 
number of modes corresponds in fact to the minimum order 
of modes to correct, as we assumed a perfect AO system 
in our study. In a real AO system, other sources of error 
like time delay, measurement noise and anisoplanatism are 
going to degrade the performance of AO, meaning that 
a larger number of modes will be required to reach the 
same ensquared energy values and image quality. We must 
especially insist on the fact that anisoplanatism is going 
to become the most important limitation in the case of 
a real instrument using both AO and 3D spectroscopy. 
Indeed, as extragalactic studies require to work at high 
galactic latitudes, the probability of finding a suitable star 
to perform wavefront sensing is equal to a few percents, 
meaning that AO, at least in its classical form, cannot be 
used because of its very low sky coverage. 

We are therefore going to describe in the next sections 
the principle and the performance of FALCON, a project 
of a multi-object 3D spectrograph with AO, which solves 
the sky coverage's problem imposed by classical AO. FAL- 
CON uses the principle of Multi-Object Adaptive Optics 
(MOAO). Such a technique allows to widen the useful field 
of view of AO systems, thus to perform the simultaneous 
3D spectroscopy of several distant galaxies in a very wide 
field of view, providing therefore a huge gain in terms of 
observing time efficiency. 



(SNR), allowing to perform 3D spectroscopy of distant 
galaxies, thus studying their kinematics. We have then per- 
formed some simulations allowing to give some preliminary 
performance of an integral field spectrograph providing an 
angular resolution of 0.25 arcsec and a spectral resolution 
R — 10000, and which would observe distant galaxies 
located at 2: = 0.9 and z — 1.5. In that case, the Ha 
emission line is rcdshiftcd in the central wavelengths of 
J and H bands (respectively 1.25 jj,m and 1.65 pm). To 
quantify the performance of such an instrument, we have 
defined a new critcrium called "^Effective signal to noise 
ratio'" (ESNR). We also assumed for those simulations 
galaxies with properties consistent with the ones observed 
in distant surveys. We have therefore shown that an ESNR 



4 THE FALCON'S AO SYSTEM: PRINCIPLE 

4.1 Increase of the AO corrected FoV: MCAO 

and tomography 

As stated in the introduction. Adaptive Optics, in its clas- 
sical form where one guide star (GS) and one deformable 
mirror (DM) conjugated to the pupil are respectively used 
to measure and correct the turbulent wavefront, suffers from 
a very low sky coverage (less than 5%), even at galactic 
latitudes of 30°. The reason of this limitation is that the 
probability to find a suitable GS in the isoplanatic patch is 
generally very low, due to the distribution of the stars in 
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our galaxy. 

Such a limitation has leaded to the development of 
Multi-Conjugate Adaptive Optics (MCAO). In MCAO, 
several deformable mirrors conjugated to the turbulent 
layers are used to correct the phase in the turbulent 
volume. Indeed, this is the repartition of the turbulence in 
the altitude which is responsible for anisoplanatism, as the 
wavefronts coming from different directions do not cross 
the same volume of the atmosphere (especially for high 
turbulent layers), and therefore do not suffer from the same 
degradations. By correcting the phase in the turbulent 
volume, and especially in the strongest layers, it is possible 
to correct the wavefront for any direction, and to have an 
uniform compensation in a field of view larger than the 
isoplanatic patch. 

Before correcting the phase with the different DMs, 
it is required to know the perturbation in each layer to 
apply to each DM the adequate c ommands. This can be 
achie v ed thanks to t omography l|Ragazzoni et al" I 1 19991 . 
I2OOOI : iTokovinin fc Via rd 2001), where the light coming 
from several off-axis GSs far outside of the isoplanatic patch 
is used to probe the 3-dimensional phase perturbations in 
the atmosphere, generally by solving an inverse problem as 
in medical imaging. However, thanks to tomography, it is 
also possible to know the integrated phase perturbation in 
the pupil for any direction in the field of view. This means 
the field where GS can be fou nd is widened, and this allows 
a higher sky coverage for AO l|Tokovinin et aLlbOOll ). 

However, the potential extension of the field of view 
with MCAO encounters some technical limitations, in par- 
ticular on the optica l system. As a n example, the ESO-MAD 
system for the VLT iMarchetti et al. (200^) or the Gemini- 
South MCAO system lEUerbroek et al.l (|20o3 ) will deliver a 
corrected field of view with a diameter of 2 arcmin, definitely 
smaller than the 10 x 10 arcmin^ field of view required for 
extragalactic studies. We are therefore going to show in the 
next paragraph a new approach for AO, where the goal is 
not to provide an uniform correction in a very wide field 
of view, but only for specific areas into it, i.e. the galaxies 
on which we want to perform 3D spectroscopy. This new 
approach is called Multi-Object Adaptive Optics (MOAO). 

4.2 Multi-Object Adaptive Optics 

First proposed by iHammer et al.l (I2OO2I '), Multi-Object 
Adaptive Optics (MOAO) is a new method of AO correc- 
tion which has essentially been since t his date the object of 
a few conference proceedings papers ijAssemat et al.l |2004| : 
iDekanv et al. 1 12004 iHammer et al.|[20o3 r We summarise its 
principles in this section. 

Let us recall here that our goal is to measure the 
internal kinematics of distant galaxies spread over the 
VLT Nasmyth field. In other words, this means that a 
high angular resolution is required only for the scientific 
targets. As a result, we propose in this section a totally 
new approach for AO: instead of correcting the whole 
10 X 10 arcmin"^ field (which anyway is impossible to do), 
we propose to correct locally only the regions of interest, i.e. 



the integral field units positioned on the scientific targets, 
and we have therefore one AO system per integral 
field unit. As a result, if we suppose that we want to 
perform the simultaneous 3D spectroscopy of 20 galaxies, 
we arrive to the concept of an instrument using 20 multiple 
AO systems, spread in the VLT Nasmyth focal plane, and 
working in parallel independently of each other. 

Although the complexity of such an instrument may 
appear insuperable at first sight, it can be greatly relaxed 
thanks to the tininess of the compensated field required for 
each independent galaxy. Inde ed, as already demo nstrated 
by the instrument GIRAFFE (|Hammer et al.lll999l '). a field 
of view of only 3x2 arcsec? is sufficient to measure the 
velocity fields of large spirals at z ^ 1, as well as galaxies 
with smaller apparent sizes observed at greater redshifts. 
This is one order of magnitude smaller than the isoplanatic 
patch in the near-infrared, so only one unique deformable 
mirror conjugated to the pupil is required to correct the 
whole galaxy. 

Moreover, we propose a concept where the AO system is 
miniaturized. Each IFU, laying in the telescope focal plane, 
will include its proper micro deformable mirror -with a 
proper tiny pupil imaging optics. The wavefront sensors also 
lay a few centimeters apart, in the focal plane ; they are lo- 
cated on the guide stars surrounding the galaxy. Translated 
into length units in the VLT focal plane, an apparent size of 
3x2 arcsec? covers a physical size of 1.8 x 1.2 mm'^ , and a 
median distance of 1 arcmin between galaxy and guide stars 
translates into ~ 35 mm. Those numbers suggest the size 
of the system we propose: the adaptive IFU, and the WFSs 
have to be integrated into suitable optomechanical devices, 
which should not be too wide (typically less than 20 mm) 
in order not to obstruct the focal plane and allow to use of 
the closest suitable guide stars for wavefront sensing. Such 
an integration of the subsystems, together with a separation 
of their functions (the DM and the WFS become separate, 
physically independent items), allow to move those devices 
into the VLT focal plane just as the positioner OzPoz does 
with the IFUs used on GIRAFFE (Hammer et al. 199^. 

In addition to offer an obvious multiplex advantage, the 
proposed structure combined with a tomographic wavefront 
sensing approach allows us to overcome the sky coverage 
problem. Several guide stars, far out the isoplanatic patch, 
are used around each IFU to sense the wavefr ont. For a 8 
meter diameter telescope, it has been shown (|Fusco et al.l 
119991 ) that 3 WFSs are sufficient to reconstruct the phase in 
the pupil for any direction in a field of view of « 2 arcmin 
diameter and for any turbulence profile; using more WFSs 
only brings a marginal improvement. 

We chose to focus our study on natural guide stars. 
We believe that using Laser Guide Stars (LGS), although 
very attractive and apparently drastic solution at first sight, 
could be technically exponentially difficult to implement 
when dealing with an increasing number of beacons (mini- 
mum of 20 beacons here), while bringing its set of inherent 
problems such as cone effect, beacon elongation, superimpo- 
sition of Rayleigh scattering and beacons, etc. Our approach 
was then to start the study on NGS, at the risk of doing the 
study again with LGS in the case where it would be demon- 
strated that NGS do not allow to fill the requirements. 

We emphasize that the architecture described above 
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corresponds to an open-loop system, since the WFS does 
not get any optical feedback from the deformable mirror 
of the adaptive IFU: this drastically differs from classical 
closed-loop AO systems. This is of course extremely chal- 
lenging in terms of AO components. Firstly, the drifts, the 
linearity, the hysteresis and the calibration of the DM must 
be kept to a very low error level in order to be able to con- 
trol it accurately enough in open-loop. Secondly the WFS 
must have a high dynamic range in order to provide reliable 
absolute open-loop measurements, together with a high sen- 
sitivity versus flux level. Such a WFS does not currently 
exist and needs to be developped, although some similar at- 
tempts h ave been done in the field of a posteriori speckle 
imagmg (|Lane et all 120031 ). Thirdly, the miniaturisation of 
the components will undoubtly demand a large technical ef- 
fort. 

On the other hand, an architecture where the optical 
system is distributed all over the field has the considerable 
advantage to simplify the optics. In particular, the field 
selection functionnality becomes here a built-in function, 
when the amount of technical difficulties and the sharp 
specifications are often underestimated in classical systems. 
Not only this critical part does not exist any more, but the 
overall optical throuput of the instrument is boosted, as, at 
any moment, no dichroic plate is needed to split the light, 
and because the number of optical surfaces is reduced to 
its strict minimum. It is important to keep in mind that 
the efforts put in this instrument are undertook to gain 
not only in spatial resolution, but also in signal-to-noise 
ratio. On this latter particular point, the instrument 
optical throughput does matter even more than the AO 
performance, which, in return, will be increased when the 
WFS optical transmission is better. 

In order to evaluate the system performance in terms 
of its dimensionning, we ran some numerical simulations. 
Their presentation and results are the scope of the following 
section. Those simulations assume that we have 3 natural 
guide stars per IFU, an a tomographic reconstructor. They 
are independent of any detailled system configuration, such 
as the miniaturisation of components or the error inherent to 
open-loop aspects. In this way, they could still apply what- 
ever the technical solution. As an example, a system where 
the Nasmyth field would be optically sliced into sub-fields, 
themselves sent to independent closed-loop mono-mirror to- 
mographic multi-analysis systems, would lead to the same 
results. 



5 THE FALCON'S AO SYSTEM: 
PERFORMANCE 

5.1 Introduction 

We show in this section the expected performance of a 
MOAO system like FALCON using atmospheric tomogra- 
phy methods on natural guide stars to reconstruct the on- 
axis wavefront from any geometry of off-axis measurements. 
We wrote a Monte-Carlo simulation code able to compute 
the long exposure PSF corrected by FALCON, for any at- 
mospheric conditions (seeing, outer scale, profile) and for 
any wavelength. This PSF was then used to compute the 



Strehl Ratio, the full width at half maximum (FWHM) and 
the ensquared energy into a square aperture of any size. 

For all the reasons explained in 13.11 we have focused 
our studies on the evolution of the ensquared energy in an 
aperture of 0.25 x 0.25 arcsec? (corresponding to 1 kpc for 
1 ^ 2 3). This evolution was studied 

• as a function of the size (number of degrees of freedom) 
and sensitivity (limiting magnitude) of the AO system 

• for two wavelengths: A = 1.25 and A = 1.65 /im 
(central wavelengths of J and H bands). As we saw previ- 
ously, such wavelengths correspond to the observation of the 
Ha emission line redshifted at z = 0.9 and z = 1.5. 

• for three galactic latitudes b ~ —90°, b ~ —60° and 
b ~ —30°, in order to study the gain brought by atmospheric 
tomography methods as a function of guide star density, and 
estimate the corresponding sky coverage, 

• for the median atmospheric conditions of the VLT site 
(Cerro Paranal, Chile). 

A particularly important aspect of this code is the sim- 
ulation of the tomographic reconstruction, this latter being 
detailed in the section [5^ 

We did not want to make technical assumptions about 
the type of the deformable mirror. For the sake of simplicity, 
we chose its infiuence functions to be Zernike polynomials. 
For the same reasons, we assume that the wavefront sensor 
measurements directly give the coefficients of the Zernike 
expansion of the measured wavefront. 



5.2 Optimal reconstruction of the on-axis phase 
from off-axis measurements 

We focus here on the reconstruction of the on-axis galaxy's 
phase in the pupil, from the off-axis wavefront sensors 
measurements. The goal is therefore to use NGS for 
wavefront sensing, which in that case are very likely to 
be located far outside of the isoplanatic angle So, and to 
determine the best commands to apply to the DM from 
off-axis measurements, in order to correct the wavefront 
coming from the galaxy's direction. 

We recall in this section the main points of tom ographic 
reconstruct ion, as developpe d by Fusco et al.l (|200l[ ). himself 
inspired bv lRagazzoni et al.l l|l999h . At a given altitude the 
phase perturbation is described by its expansion on a modal 
basis, defined over a disk wider than the pupil a n d cov ering 
any beam imprint of the field l|Ragazzoni et al] l|l999l ) call 
metapupil this area). Assuming a finite number of turbulent 
layers, the phase in the whole volume is reduced to one vec- 
tor, and the game is to express the phase expansion within 
the pupil, for a given direction, as a matrix product with this 
vector. Once this relation has been written, the inversion of 
this linear relation has various solutions, depending on the 
approach followed by the authors (least-square, maximum 
likelihood, etc). 

Let us discretise the turbulence profile so that it can be 
modeled as a finite number Nt of turbulent layers. Assuming 
near-field approximation, the resulting phase at a position r 
in the pupil for a sky direction cti can therefore be written 
as the sums of the phase ipj produced by each turbulent j'"" 
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Figure 11. The Falcon concept combines MOAO and tomography: several IFUs (rectangles) are spread over the VLT focal plane. 
Each IFU is used to perform the 3D spectroscopy of one galaxy. Three miniaturised WFSs (circles) measure the wavefront coming from 

3 ofT-axis NGS located around the galaxy. Their combined signals are used, through a tomographic reconstructor, to command the 
miniaturised DM into the IFU and correct the wavefront. It is therefore possible to recover the velocity field of the galaxy with a better 
spatial sampling. 



layer at the altitude hj : 

*(r,ai) = ^ipj{r + hjoti) 



(6) 



If our system of coordinates is centered on the galaxy 
(on-axis, oci = 0), then the phase perturbation from the 
galaxy will be written as: 

$G(r,0) = 53^,(r) (7) 

3 = 1 

Each phase $(r,ai) can be expanded as a sum of 
Zernike polynomials where we omit the piston term Zi(r): 



(8) 



We assume now that we sense the wavefronts coming 
from a number Nqs of stars around the galaxy, i.e. we have 
Nas WFSs, and that we use the combined measurements of 
all those WFSs to compute the on-axis phase coming from 
the galaxy. Moreover wc assume that the WFSs measure- 
ment are the coefficients 2 to kmax of the Zernike expansion 
of the phase, plus some noise. We will now call the measured 
phase the following quantity: 



$"'(r,ai)= aa,,kZk{r) + ni{r) 



(9) 



where ni(r) is the measurement noise for the WFS look- 
ing at the direction cxi, i.e. the propagated noise on the 



Zernike polynomials in the wavefront reconstruction process. 
We can therefore write: 



(10) 



The m^k coefficients are stored in a rii vector, which 

is a random gaussian vector whose statistics is given by its 
covariance matrix C„,i, dependent of the type of WFS as 
well as GS magnitude and readout noise. 

We call now ^ai,kmaa: the vector storing the coordi- 
nates of $(r,ai) up to the Zernike polynomial kmax'- 



\ 0,ai,kmax I 



(11) 



This allows us to define the vector storing the co- 
efficients of the phase 3>'"(r,ai): 



aai,3 + Wi,3 



*ai,fe^„„+ni (12) 



Let us consider now the phase perturbation 'Pj{Pj) in 
the i*'^ turbulent layer. This latter can be expressed also as 
a sum of Zernike polynomials. We call ipj the vector storing 
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the coefScients of this phase up to the index Nj, and we 
have: 



aj,2 
aj,3 



(13) 



If we concatenate now all the vectors tpj, j — 
l,2,...,Nt, we have the vector (p storing the coefficients 
of the phase in the volume: 



[Ragazzoni et al.' 



¥=2 



19991): IFusco et all 



(14) 



iFem enia & Dcvancv (2003) have shown that there is a 
linear relation, that we will not detail here, between 
^ai.femca, and ip. It takes the form: 



(15) 



where the matrix is called a "modal projection" 

matrix, performing the sum of the contributions of each 
wavefront Lpj{pj) on the telescope pupil for a given direc- 
tion cti. Therefore, if we concatenate all the vectors, 
i = 1,2, ... , Nqs into the single vector we can write 

ijFusco et al.ll200ll ): 



(16) 



where n is the concatenation of all the measurement 
noises rii, i = 1, 2, . . . , Nqs, and M^* is the generalisation 
to several directions of the matrix Ma^ ■ 

Let us consider now the on-axis galaxy's phase ^cir, 0). 
This latter can also be written as a sum of Zernike polyno- 
mials: 



$g(i-, 0) = ^ aG,kZkir) 



(17) 



Let $G be the vector storing all the aG,fc coefficients up 
to oo, and ^Cfemaaj be the vector storing these coefficients 
up to the Zernike polynomial kmax- This last vector is also 
linked to the phase in the volume ip by the following relation: 



(18) 



where Mg ' is the matrix summing the contributions of 
each wavefront 'Pj{Pj) on the telescope pupil for the on-axis 
direction ai = 0. Our goal is to find the best estimation 
^G,femaa!5 whlch wlU be the vector storing the commands 
to be applied to the DM, provided that its infiuences func- 
tions are also the Zernike polynomials 2 to kmax- Therefore, 
the correction phase $G(r,0) provided by the DM has the 
following expression: 



$G(r,0)= ^ aG,kZkir) 



and we seek a relation in the form: 



(19) 



(20) 



where W is the reconstruction matrix linking the com- 
mands ^G,k^a.a! to the off-axis measurements Let us 
define a minimum mean-square error (MMSE) for the re- 
construction. In our case we want to minimize the residual 
variance of the wavefront ^es ■ 



^^air) - $G(r')] dr 



W 



(21) 
(22) 



1 

^ J 

= £ ((aG.fe-aG..)'>+ J2 fe> (23) 

fc = 2 fc = fcmax+l 

This residual variance is therefore the sum of two terms: 

• X^fc=2^ {{cLG,k —ac.k)^), the reconstruction error due 
to anisoplanatism and noise propagation 

• X^fc^fc (flG.fc), the uncorrected variance due to the 
finite numbers of Zernike polynomials used to correct the 
wavefront. 

We therefore want to minimize 

X]fc=2^ (('^G.fc — ^G.fe)^)- The optimal recon struction 
matrix satisfying this condition can be written as (|Wallnej 
1983): 



Wopt = ( * 



(24) 



i.e. this matrix is the product of two matrices: the co- 
variance matrix of the unknowns and of the measurements, 
and the inverse of the covariance matrix of the measure- 
ments. This expre ssion is equivalent to the one obtained by 
IFusco et al.l (I2001I ) with a maximum a posteriori approach. 



Using the equations p6[) and (|18p . and as the noise is 
uncorrelated from the turbulent phase, the equation (|24|l 
can be rewritten as: 



opt 



Cb 



(25) 

We recognize in this expression the product of two 
terms: 

• the matrix {M^^f \mS^C^ (M^')^ + Cb] \ 
which gives the best tomographic estimation of the phase 
in the volume tp from off-axis measurements. We find here 
the same expression than in the equation (18) of lFusco et al.l 
(2001) 

• the matrix JV/q *: as explained before, it sums the 
contributions of each wavefront fj[Pj) on the telescope 
pupil for the on-ajcis direction a; = 0. 

An important point here is the presence of th e matri- 
ces and Cb, introduced bv IFusco et al.l l|200ll ). which 
are the generalisation for several layers and for several GSs 
of the classical turbulence and noise covariance matrices. 
Thanks to the information contained in those matrices, it is 
possible to regularize the inversion and increase the field of 
view where off-axis NGS can be picked off to perform wave- 
front sensing, as well as to use fainter NGS than in classical 
least-square methods. However those matrices require some 
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a-priori knowledge on the noise measurement, and on the 
turbulence profile. 

On a practical point of view, equation (I24|l shows 
that "half of the optimal matrix can be measured in-situ: 
the covariance matrix of the measurements is something 
that can be obtained from the experiment itself. Only 
the covariance matrix between real, actual phase and 
measurements has to be computed. 

It must be noticed that we used Zernike polynomials 
in our simulation because they are easy to use, useful 
mathematical tools. However the current study can be 
applied to any other modal basis. Further work about 
the optimal reconstruction matrix in the case of real AO 
components (Shack-Hartmann or pyramid WFS, segmented 

or continous facesheet DMs) has already been started 

l|Assemat.,2004i '). 



5.3 Detained presentation of the simulation code 

The elements of this code are the following: 

• a mono or multi-layer turbulent atmosphere. The num- 
ber of layers, their altitude, the outer scale £o as well as 
the strength of the turbulence in each layer (given by a lo- 
cal ro) are adjustable. The atmospheric p hase screens ar e 
simulated using Fourier filtering methods l|ShaklanllT989l ). 
and take into accounts the effects of the outer scale Cq by 
introducing a proper Von-Karman spectrum 

• any number of NGS, whose position in the field as well 
as magnitudes are adjustable 

• any number of WFSs positioned on the off-axis NGSs. 
The WFS is supposed to give the Zernike expansion of the 
phase. Some noise is added to the true Zernike coefficients. 
The noise is related to the guide star magnitude using the 
propagation properties in the reconstruction of Zernike poly- 
nomials. It is then possible to have Nqs measurement vec- 
tors , which are concatenated in the vector 

• the computation of the reconstruction matrix Wopt 
from the GS geometry and magnitudes, the turbulent pro- 
file and the number kmax of measured and corrected Zernike 
polynomials 

• the computation of the vector ^G,kma.^ = Wopt 
storing the coefficients of the reconstructed on-axis phase on 
the Zernike polynomials 2 to kmax 

• a DM with {kmax — 1) actuators, whose influence func- 
tions are also the Zernike polynomials 2 to kmax , and whose 
commands are stored in the vector ^G,fem,<iaii allowing to 
compute the corrected phase I>g(t',0) = J^'k^T o,G,kZk{r) 

• the computation of the on-axis residual phase 
^res{r, 0) = ^G{r) - ^G{r) 

• the computation of the short exposure AO corrected 
PSFs at different imaging wavelengths 

• the computation of the long exposure AO corrected 
PSFs by averaging the AO short exposure corrected PSFs. 

We considered in all those simulations a 8 meter tele- 
scope (VLT case), whose pupil was simulated on a discrete 
grid of 128 x 128 pixels. Assuming Nyquist sampling in the 
focal plane (A/2D), this leads to simulated PSFs covering 
respectively a fleld of 4.12 x 4.12 arcse(? in J band and 



5.44 X 5.44 arcse(? in H band. Each long exposure AO 
corrected PSF was the average of 100 independent short 
exposure PSFs, as we assumed an open-loop system and 
that we did not consider any temporal error. 



5.4 Simulation conditions 

We are giving in this paragraph more details about the con- 
ditions of the simulations we performed. 

5.4.1 Atmospheric conditions 

Let us first deal with atmospheric conditions. We used data 
from ESO AO's department, who provided us with some 
statistics about the seeing &1/2, the isoplanatic angle Oq 
and the coherence time tq. We therefore considered some 
median atmospheric conditions, leading us to the following 
quantities (all data given for zenith and for a wavelength of 
0.5 /im): 

• a median seeing 6*1/2 = 0.81 arcsec. This latter is linked 
to the Fried parameter ro by the relation 6*1/2 = 0.976 A/ro. 
We therefore find a median Fried parameter ro = 12.42 cm 
at 0.5 iim for the whole turbulence profile 

• a median isoplanatic angle do — 2.42 arcsec 

• a median coherence time ro = 3.04 ms 

We therefore used the median seeing and isoplanatic 
angle to define our turbulence profile, and find that a profile 
made of 3 turbulent layers located at altitudes of (ground 
layer), 1 and 10 km, and with respectively 20%, 65% and 
15% of the whole turbulence allowed to reproduce these 
atmospheric conditions. 

Another important issue is the outer scale Co, which 
has a direct infiuence on the variance of low orders of 
the turbulence a nd the FWHM of uncorrected images 
(!Tokovininll2002l ). iMartin" et al. I (I2OOOI) give some statistics 
for this parameter at the Cerro Paranal, and we adopted 
their median value, i.e. £.0 ~ 24 m. 



5.4.2 AO system parameters 

We focus now on the components of the AO system. We 
considered for each simulation case a tomographic system, 
with 3 off-axis wavefront sensors and one deformable mirror. 
As said before, we assumed a correction degree ranging from 
to 120 Zernike polynomials. We assumed we use wavefront 
slope sensors (Shack-Hartmann or pyramid), leading to a 
propagated noise v ariance on Zernike polyno mials following 
a law in (n + 1)"^ (|Rigaut fc GendronllToQ^ ). 

5.4.3 AO system limiting magnitude 

The result of our study depends on the WFS noise level, 
and this latter is a function of the photon flux. We assumed 
to be in a regime dominated by photon noise. In that case, 
the noise variance is invers ely proportional to the nu mber of 
photoelectrons per frame l|Rigaut fc GendronI Il993 ). N~i^ . 
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Hence, the noise behavior of the sensor will entierely be de- 
fined versus flux and for any Zernike mode, when the con- 
stant of proportionnality of the law is given. 

The way we propose to fix this constant of proportion- 
nality is to define the flux level that produces a particular 
noise variance. In order to make this flux level number 
physically meaningful, we will choose for our noise variance 
value, the one which corresponds to the limiting magnitude 
of the system. Doing this, setting a limiting magnitude 
implies a particular noise variance, and consequently deflnes 
the complete noise behavior for any flux level. The usual 
way of doing it is to derive the WFS performance (and the 
limiting magnitude) from the list of the WFS characteristics 
: detector noise, optical throughput, spectral bandpass, 
number of lenslets, etc. We think it is not reasonnable 
to undertake such a study : not only because this would 
always be possibly technically quest ionnable, out of the 
scope of the article, but also because a study based on 
current existing designs might be obsolete in a few years. 

Instead, our approach is to consider the WFS limiting 
magnitude as the input parameter to our study, not as 
an output number resulting from a complete modelisation 
of the wavefront sensor. Here, it should be seen as a 
specification number, when it's a matter of designing the 
wavefront sensor. The question is then to decide of a 
value for the noise variance that characterizes the limiting 
magnitude. The choice is just a matter of definition. In 
order to give meanigful numbers, we took a criterion 
that agrees with actual AO-users' experience. Practically 
speaking, the criteria is reached when the user might be 
discouraged trying higher magnitudes, considering that the 
correction only brings marginal improvement. From our 
experience with the NAOS Adaptive Optics system, behind 
this fuzzy rule of thumb is hidden a clear-cut number : we 
choose a criterion of a noise level value of 250 rd^ at 0.5 
microns. This number really fits with what is actually called 
"limiting magnitude" by the AO community. It cor responds 
to a magnitude oi R = 17 (|Fusco et al.l |2004 ) on the 
NAOS system (whose detector on the WFS has a read-out 
noise equal to 4 e- rms per pixel). At this magnitude, the 
NAOS system uses 7x7 subapertures and controls 42 mir- 
ror modes, only 10-15 of them being efficiently compensated. 

Now, the choice of the value for the limiting magnitude 
has to be done. We cannot arbitrarily increase it, and in 
particular we can set an upper limit, translating into lim- 
iting magnitude the classical laws of photon-noise limited 
wavefront sensing, that are widely used in the litterature. 
Following iRigaut fc GendronI (|l992l ) , when the subapertures 
are larger than ro, the noise variance can be written as: 



with no the number of photoelectrons per surface area. 
Using that formula, if we consider a spectral interval of 
0.4 fim, a sampling frequency of 100 Hz, and a global 
optical throughput of 0.5, then the value — 250 rd^ is 
attained when no — l.QZphe" /m? / frame, i.e. at a limiting 
magnitude of _R = 20.0. This number must be taken as 
the very upper limit, since this theoretical formula only 
assumes quantum-noise limitation, excluding any other 



instrumental effect. 

We know that in a near future, progress in the fast 
readout detectors may lead to m uch better results than 
NAOS. As an example, L3CCD (iDaigle et all 120041 ') are 
extremely good candidates, that could bring the noise 
level below 1 e~ rms per pixel, i.e. increase the limiting 
magnitude of -1-1.5 compared to NAOS. However, we have 
decided to be rather conservative, and we will study in this 
article two cases of limiting magnitudes, Rmax ~ 16 and 
Rmax = 17, that we think to be realistic. In that case, one 
may then wonder how our results shown in section 15.71 can 
exhibit corrections up to 120 Zernike polynomials, when 
we noticed that a sensor working at its limiting magnitude 
-considering our definition- can only efficiently compensate 
for 10 to 15 Zernike modes (out of 42 effectively handled, 
and partially corrected) . The answer is that the magnitudes 
of the 3 guide stars are not simultaneously equal to Rmax'- 
as shown in the figure [121 the star surface density increases 
as a function of the R magnitude by a factor 2 every 
1.5 magnitudes for magnitudes R > 13. It means that 
the median magnitude of a field of stars brighter than 
R = Rmax is equal to Rmax — 1.50. Then as the signal of 
the 3 wavefront sensors are combined together through an 
optimal reconstructor where both statistics of measurement 
noise and anisplanatism are taken into account, thus 
increasing the final SNR compared to the one that would 
be given by a single WFS, this still allows us to partially 
compensate for high-order modes. And as it was explained 
in section 13.21 such a correction is particularly effective 
at increasing the ensquared energy of the AO corrected PSF. 



5.4-4 Star density 

Let us consider then the three off-axis NGS used to per- 
form off-axis wavefront sensing and on-axis wavefront re- 
construction. The angular distance of the first, second and 
third closest NGS to the central object, is directly linked to 
the galactic latitude, a s the surface density of st ars decrease 
with galactic latitude (|Bahcall fc Soneiralll980l ). and has a 
direct impact on the performance of the tomographic AO 
system. In order to have some estimation of the sky cover- 
age of the falcon's AO system, we therefore decided to 
perform our simulations on three stellar fields, located at 
the following galactic coordinates: 

(i) / = 280.2°, b= -88.5° 

(ii) I = 0.3°, fo = -61.4° 

(iii) / = 44.5°, b = -30.9° 

We us ed statistics given by the Modele de BesanQon 
dRobin e"t~ al. 2003) to simulate star flelds for these three 
galactic latitudes. This requires to know the surface density 
of stars, which are shown on the figure 1121 and the table 
[2] shows the integrated number of stars per square degree 
for the limiting magnitudes Rmax ~ 16, Rmax = 17 and 
Rmax = 18. Then, for each field, we choose 100 random 
positions. For each position, we chose the three closest 
NGS with a magnitude R ^ Rmax to perform wavefront 
sensing, and computed the reconstruction matrix given 
by equation ((25} for each maximal degree of correction. 
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b 


N{Rma^ = 16) 


N{Rma^ = 17) 


N{Rmax = 18) 


-30° 


1271 


2048 


3096 


-60° 


400 


639 


946 


-90° 


246 


361 


522 



Table 2. Number of stars per square degree as a function of the galactic latitude b and the magnitude Rmax 




R mag 

Figure 12. Integrated star counts per square degree in R band as a function of the R magnitude given by Modele de Besangon for three 
galactic latitudes: b -90° (full hne), b -60° (dashed line) and b Ri -30° (dotted line). 



We were then able to reconstruct the on-axis phase from 
off-axis measurements, compute the on-axis residual phase, 
and the long exposure AO corrected PSF for different 
correction orders in J and H band. 

The quality of the on-axis wavefront reconstruction is 
directly linked to the geometry of the closest NGS used to 
perform wavefront sensing: anisoplanatism will be increased 
when the NGS will be far away of the scientific target. More- 
over, as we choose the three closest NGS, it is very likely that 
those latest will be faint, meaning an important measure- 
ment noise. We are therefore giving in the next paragraph 
some analytical formulae for the distances of the first, second 
and third closest NGS. 



5.5 Distances of the closest NGS 

Let us call a* the surface density of stars, i.e. the number 
of stars per angular surface on the sky. As an exemple, the 
previous paragraph has shown that we have cr, — 246 stars 
per square degree for Rmax = 16 and b ~ —90°, thus a, — 
0.07 stars per square arcmin. Therefore, if we call A/'(r) the 
average number of stars in a circle of radius r (in arcmin), 
we'll have Af{r) — ivatr^. As the spatial distribution of stars 
on the sky follows an uniform probablity law, the number of 
stars in an area of 7rCT*r^ arcmin^ follows a Poisson law. So, 
the probability to find k stars in an area of yrr^ arcmin^ is 
equal to: 



nC = fc) = -^exp( 



(27) 



The probability that there is then one star at a distance 
greater than r is also equal to the probability that there 



is no star in an area equal to nDr^ arcmin^, i.e. P(^ = 
0) = exp —Af{r) = exp (— 7r(7«r'^) . Then, the probability that 
there is one star at a distance _Di ^ r is equal to 



Pi(r-) = P(Di ^ r) = 1 - P(e : 
= l-exp(-AA(r)) 
= 1 — exp (— 7r(T,r^) 



0) 



(28) 



and we recognize the cumulative distribution function 
(cdf ) Pi (r) of the random variable Di , this latter being the 
distance of the first closest NGS. Then, the probability of 
having one star at a distance between r and r + Ar will be 
equal to the derivative of the cdf of Di, i.e. its probability 
density function (pdf) pi{r), which is equal to: 

pi(r) = '^^^^'^^ 



= 27r(j,r exp (— A/'(r)) 
= 27r(Jtr exp (— TTtT.r^) 



(29) 



As an example, the figure [TH] shows the cdf and the 
pdf of Di for a galactic latitude b ~ —90° and a limiting 
magnitude Rmax = 16. Knowing it, it is possible to know 
two interesting quantites. First the average value of Di: 



r pi{r)dr 



'0 

1 



2 / 

r exp I 



(30) 



and the most likely value of Di, that we'll note Di^mi, 
given by the maximum of the pdf Pi{r). The determination 
of this maximum requires to compute the derivative of pi (r ) , 
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Figure 13. Cumulative distribution function (left) and probability density function (right) of the distance of the first closest NGS 
computed from the star density at the galactic pole {b —90° , Rmax = 16) provided by Modele de Besangon. 



whose expression is; 

—pi{r) = 27r(T, exp (— 7ro"*r^) (l — 27rcrt 
dr \ I \ 

and which is equal to for 



(31) 



(32) 



^/2^ 

Let us focus now on D2, the distance of the second 
closest NGS. Doing the same assumptions than before, the 
probability to find two stars at a distance D2 ^5 r is equal 
to: 

P2{r) = P{D2^r) 

= i_[p(e = o) + p(^ = i)] 

= l-exp(-AA(r))(l+AA(r)) 
— 1 — exp (— TTfT.r^) (1 + 7r(T*r^) (33) 
The pdf of D2 has then the following expression: 
P2{r) = 27r(j»rA/'(r-) exp (— A/'(r)) 



o 2 2 3 I 2\ 
= ZTT (7,r exp (— 7ra,r j 

and the average value of {-D2) can be computed: 

f-oa 

{D2) ^ r p2{r)dr 

Jo 

= 2na, / r^A/'(r)exp(-A/'(r)) 
Jo 

3 



The derivative of P2{r) is equal to: 

■^P2 (r) = 27ra. AA(r) exp -7V(r) [3 - 2f^(r)] 
dr 

= 2TT^a^r^ exp — 7r(T*r^ [3 — 27ra,r^] 
and is equal to for 

D2,ml - 



(34) 



(35) 
(36) 



(37) 
(38) 



We can now focus on D3 , the distance of the third clos- 
est NGS. The probability that there are three stars at a 
distance Dj, ^ r is equal to: 

PaW = P(I>3 =S r) = 1 - [P(5 = 0) + P(e = 1) + P(C = 2] 



-exp(-AA(r)) l+AA(r) + 



2 2 4 

^ 2 1" O"*'' 

1 — exp (— 7r(T,r ) ( 1 + 7ra»r + 



The pdf has then the following expression: 



(39) 



p-i{r) = 



335 
= TT a^r exp 



rjV^{r)exp{-jV{r)) 

2\ 



(40) 



and can be used to compute the average distance {D3): 

/oo 
r pz{r)dr 

r^7V^(r)exp(-7V(r)) 



Jo 
7ra» 



15 



16v^ 

The first derivative of pz(r) is equal to 



(41) 



d 



= Tva*Af'\r) exp -A/'(r) [3 - 2Af{r)] 



and is equal to for 

Ds,mi = 



2Tva, 



(42) 



(43) 



The table [3] summarises the above results, whereas 
the table |4] gives some values of distances for the three 
galactic latitudes we studied and for the two cases of 
limiting magnitude. We can notice there that the distances 
of the closest NGS fastly increase with the galactic latitude 
b. As an example, for a limiting magnitude Rmax = 16, 
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(D) 




pi 


27ro'tr exp (— 7rcr,r^) 


1 






27r^(T^r"^ exp (— 7r(T«r^) 


3 

4^ct7 


/ 3 
V 2irCT, 


P3 


TT'^crJr^ exp (— Trcftr-^) 


15 


\/ 2ir(T» 



Table 3. Summary table for the probability density function, the average and the most likely value of the distance of the closest first, 
second and third NGS used for tomography as a function of the star surface density cr. 



b 


Rmax 








{D2) 


D2,ml 






-30° 


16 


0.35 


40 


50 


70 


76 


90 


94 


-30° 


17 


0.57 


32 


40 


55 


60 


71 


75 


-60° 


16 


0.11 


72 


90 


124 


135 


160 


169 


-60° 


17 


0.17 


57 


71 


98 


107 


127 


133 


-90° 


16 


0.07 


92 


115 


158 


172 


205 


215 


-90° 


17 


0.10 


76 


95 


130 


142 


169 


177 



Table 4. Average and most likely distance of the closest first, second and third NGS as a function of the galactic latitude and the 
limiting magnitude. Prom left to right: the galactic latitude b, the limiting magnitude Rmax, the number of stars per square arcmin ct*, 
and the average and most likely value of the distances of the off-axis NGS (in arcsec). From star densities given by Modele de Besangon. 



the distance Di,™; goes from 50 arcsec for b — —30° to 
115 arcsec at the galactic pole. The difference is more 
important for the third NGS: we have Di.mi ~ 94 arcsec 
at 6 = —30° (which remains usable thank to tomography), 
but at the galactic pole, the third NGS is very likely to be 
at a distance of more than 3 arcmin ! The gain brought by 
a higher sensitivity becomes then obvious: if fainter NGS 
(i.e. Rmax = 17) are used, these distances dramatically 
decrease. At a galactic latitude b = —30°, the distances 
of the closest first and third NGS are respectively 10 and 
20 arcsec smaller than for Rmax = 16, and for greater 
galactic latitudes (|fe| ^ 60°), these distances are 20 and 
35 arcsec smaller. Such differences are not negligible for 
AO and particularly for anisoplanatism's correction, and 
show that the improvement of WFS's sensitivities is really 
mandatory for future systems. 



5.6 Sky coverage: definition of a criterium 

One of the key issues of AO systems is their sky coverage, 
i.e. the probability to reach a certain performance thanks to 
the correction provided by AO. For the studies of new AO 
systems, it is very useful to be able to know the variability 
law of the sky coverage as a function of observing condi- 
tions, like the galactic latitude or the limiting magnitude of 
off-axis guide stars. 

This implies therefore to set the performance we want 
to reach. As we have seen in the section [3] ensquared en- 
ergies of 30% in J band and 40% in H band could allow to 
reach a sufficient SNR to perform the dynamical studies of 
galaxies located at z ^ 0.9 with an angular resolution of 
0.25 arcsec and a spectral resolution R — 10000, and a min- 
imum exposure time of 3 hours. Such values of absolute 
ensquared energy can therefore be used as a basis to design 



the instrument, and we show in the next sections the results 
we obtained for this specification. Definitive values of the 
absolute ensquared energy would however require to made 
the same simulations for a broader sample of galaxies, which 
currently are not possible because of the current lack of 3D 
spectroscopic data for distant galaxies. 

5.7 Results 

The figures [14] and [15] show the median performance (over 
100 NGS triplets) of the FALCON's AO system, i.e. the 
evolution of the ensquared energy in an 0.25 x 0.25 arcsec^ 
square aperture and the evolution of the FWHM as a 
function of the number of corrected Zernike polynomials. 
The figures on the left assume a WFS with a limiting 
magnitude Rmax = 16, the ones on the right a WFS 
with the same noise variance than the NAOS visible WFS 
{Rmax = 17). On each figure, 3 plots have been reported: 
they correspond to galactic latitudes of 30, 60 and 90 
degrees, and therefore show the performance for 50% of the 
simulation cases, corresponding to a sky coverage of 50%. 
The specification, i.e. an absolute ensquared energy of 30% 
or 40%, or a FWHM of 0.25 arcsec is also drawn on these 
figures. 

One should notice the global shape of these curves. 
For the ensquared energy curves, we can see that they still 
significantly increase, even at a high number of corrected 
modes. However, the FWHM curves quickly saturate 
to values which are always less than our specification 
{FWHM = 0.25 arcsec). In fact, as shown in the section 
13.21 high order modes need to be corrected in order to gather 
the light in the inner regions of the PSF and to improve 
the ensquared energy. Therefore, to reach our required 
specifications for the ensquared energy, a non negligible 
number of modes must be corrected. For a WFS with a 



The FALCON concept 23 




Figure 14. Evolution of the median ensquared energy in a 0.25 X 0.25 arcsec^ square aperture for each galactic latitude b and each 
limiting magnitude Rmax in the J and H band, as a function of the number of corrected Zernike polynomials. The full line corresponds 
to the median performance at the galactic pole, the dashed one to the median performance at 6 —60°, and the dash-dotted one at 
b fa —30°. The dashed horizontal line shows the specification, i.e. 30% in J band and 40% in H band. 



moderate quality, at least 100 (for J band) and 120 (in H 
band) Zernike polynomials must be corrected so that a sky 
coverage of 50% can be reached for any galactic latitude 
b. With a more sensitive WFS {Rmax = 17), allowing to 
use fainter and closer NGS, the correction order should 
be only 65 to 85 Zernike polynomials to reach the same 
sky coverage value: the better the WFS, the lower should 
be the correction order, as there is a transfer of the phase 
error between noise and anisoplanatism. The consequence 
is therefore that a rather large number of actuators would 
be required, even with the moderate angular resolution 



required. 

These curves show the median values over 100 NGS 
triplets of the expected performance expected thanks to the 
falcon's AO system. This means that the performance 
will be better than shown for 50% cases, and worse for the 
other 50%. In other words, it means that the performance 
shown here stand for a sky coverage of 50%. 

Immediately, a question arises, that is to know whether 
50% sky coverage is a sufficient number or not. This seems to 
be a critical issue, as demonstrated on the figure [16] where 
the ensquared energy has been plotted as a function of the 
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Figure 15. Evolution of the median full width at half maximum (in arcsec) for each galactic galactic latitude b and each limiting 
magnitude Umax in the J and H band, as a function of the number of corrected Zernike polynomials. The full Une corresponds to the 
median performance at the galactic pole, the dashed one to the median performance at b —60°, and the dash-dotted one atb^ —30°. 
The dashed horizontal line shows the specification, i.e. FWHM = 0.25 arcsec. 



number of corrected Zernike polynomials and for sky cover- 
ages equal to 10%, 50% and 90%. We see that for a spec- 
ification of 40%, the number of corrected modes should bo 
respectively equal to 60, 75, or more than 120 ! This way of 
thinking leads therefore to the conclusion that the design of 
the system is strongly dependent of the sky coverage. How- 
ever, this is only partially true, because a sky coverage of 
50% docs not mean that only half of the objects can be ob- 
served whereas the others cannot: it means that 50% of the 
objects are observable (some of them with an even better 
quality than expected), and most of the other 50% are ob- 
servable too, but with a moderate image quality. This latter 



could however be sufficient, as it will be in any case better 
than the uncompensated seeing-limited image. 



As an example, when assuming a median corrected en- 
squaxed energy equal to 40%, 75 Zernike polynomials have 
to be corrected. With such a correction, 95% of the galaxies 
can be observed with a performance ranging between 36% 
and 43%. In other words, the dispersion of the ensquared 
energy provided by AO correction is very low, making the 
50% sky coverage a very acceptable performance. 
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Maximum index of the corrected Zernike Polynomial 

Figure 16. Evolution of the ensquared energy in H band for the field located at b —60° as a function of the corrected number of 
Zernike polynomials, and for a limiting magnitude Umax = 17. The "A" line shows the performance reached for 10% of the simulation 
cases, the "B" one the performance reached for 50% of the cases, and the "C" one the performance reached for 90% of the cases. The 
dashed horizontal line shows the specification, i.e. an ensquared energy of 40%. 



5.8 Conclusion 

Thanks to numerical simulations, we have been able to show 
in this section the expected performance of a MOAO system 
using atmospheric tomography methods to reconstruct the 
on-axis wavefront from ofT-axis measurements performed on 
NGS, in order to provide simulaneously the 3D spectroscopy 
of distant galaxies in a 10 x 10 arcmin^ field of view. For 
any galactic latitude, our results show that an AO system 
using three very sensitive WFS (allowing to measure the 
wavefront on NGS with a magnitude R ^ 17) and one 
deformable mirror conjugated to the pupil per scientific 
target would allow to reach a spatial resolution better than 
0.25 arcaec in J and H bands for 95% of the objects, and 
to achieve an ensquared energy better than 30% in J band 
and 40% in H band for 50% of the objects. Greater sky 
coverage (up to 95%) seem to be achievable for slighty less 
values of ensquared energy. As explained before, this latter 
has not been definitely fixed, especially for galaxies with 
z ^ 1.5, as further dynamical and photometric data are 
required to constrain it. However this requires to correct 
a non negligible number of modes: deformable mirrors 
and wavefront sensors with at least 10 x 10 actuators and 
subapertures are required to achieve such performance. 

From a technological point of view, such a system is 
very demanding. In particular, the miniaturisation of AO 
components in order to avoid field obstruction is critical, 
and requires further developments. Some of them have al- 
ready started: microdeformable mirrors (|Zamkotsian et al] 
l200a) would be very good candidates to be used as 
wavefront correctors in each integral field unit. But the 
miniaturised and high-dynamic WFS architecture used in 
the falcon's AO system remains to be defined. 



Therefore, if such an instrumental concept is imple- 
mented on the VLT, it will then provide a huge gain in 
terms of observing time efficiency for the dynamical study 
of high redshift galaxies. As an example, recent obser- 
vations made by the AO-ass isted spectrograph SINFONI 
(|Forster Schreiber et al.ll2006l ) have allowed to retrieve the 
rotation curves of 14 z « 2 galaxies, with exposure times be- 
tween 2 and 8 hours. However, as SINFONI was used in its 
NGS mode, only two galaxies were observed with a spatial 
resolution better than atmospheric seeing, the other ones 
being observed with a spatial resolution of 0.5 arcsec. In 
other words we can say that SINFONI in NGS mode as a 
sky coverage of ~ 14%. Moreover the spectral resolution for 
these observations was between R « 1900 and R ~ 4500. 
These numbers (spatial and spectral resolution, sky cover- 
age) have to be compared with the ones shown in the sec- 
tions |3I3] and [S771 where we can see that the performance of 
an instrument like FALCON are far beyond the capabilities 
of SINFONI in its NGS mode. The performances of SIN- 
FONI in its LGS mode will probably be better in terms of 
spatial resolution and ensquared energy. However, the fact 
that very long exposures (see section I3.4.2|l are mandatory 
to probe the kinematics of distant sources with an 8 meter 
telescope provides an enormeous advantage to a multiplex 
instrument such as FALCON. 



6 DISCUSSION 

6.1 Comparison witli GLAO 

The previous section has shown that thanks to atmospheric 
tomography, the correction provided by a MOAO system 
like FALCON allows a huge gain in terms of ensquared 
energy improvement as well as FWHM reduction, thus 
allowing to perform the 3D spectroscopy of distant galaxies. 
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However the use of tomography is quite complex, in par- 
ticular because it requires to have a real-time knowledge of 
the turbulent profile for the computation of the turbulent 
profile. 

If we look at the turbulence profile definined at the 
section [531 we can see that 85% of the turbulence is located 
between altitudes of and 1 kilometers, so very close to 
the ground. Ground-Layer Adaptive Optics (GLAO) could 
therefore be used to impro ve the ensqua r ed energy as we ll 
as the angular resolution (|Rigautl [2003 : iTokovininI [20o3 ) . 
GLAO consists in averaging the measurements from several 
off-axis guide stars, and is normally very efficient when most 
of the turbulence is located at or close to the telescope pupil. 

The key issue here is to see if GLAO provides a cor- 
rection as effective as the one provided by tomography, as 
this would decrease the AO system's complexity. To answer 
to this question, we did again the same simulations than 
the one described in the section [5] But this time, instead of 
using the optimal tomographic reconstruction matrix Wopt 
to compute the on-axis correction phase, we simply did an 
average of the off-axis measurements coming from the three 
closest GS of the science target. 

The figures [17] and [18] show therefore the median per- 
formance (over 100 NGS triplets) of the evolution of the 
ensquared energy and the FWIfM for the three cases of 
galactic latitudes and the two cases of limiting magnitudes 
already studied, but this time for GLAO correction. For the 
FWHM, we observe also a saturation of the FWHM, but we 
can see that even the use of stars with magnitudes R ^ 17 
does not allow to reach the required angular resolution for 
all the galactic latitudes, as the median angular resolution 
reached at the galactic pole is equal to 0.30 arcsec. The re- 
sults are worst for the ensquared energy curves, which show 
that it is never possible to reach the required ensquared en- 
ergy values: a GLAO system has a 0% sky coverage for 
the FALCON specification. 

Moreover, a careful examination of the ensquared en- 
ergy curves show that their behaviour is quite different from 
the ones shown on the figures [14]: contrary to these latter, 
the curves start to increase, reach an upper limit, and then 
systematically decrease. This result has to be related with 
the study of the AO corrected PSF made at the section \3l2\ 
Indeed, we have shown there that high order modes must 
be corrected to bring the energy into the inner parts of the 
PSF. In GLAO, as an averaging of off-axis measurements is 
done, there is an important residual anisoplanatism which 
degrades high-order correction, thus explaining the decrease 
of the ensquared energy for an important number of cor- 
rected modes. This effect is moreover going to be degraded 
by the noise propagation because of the faint stars used for 
wavefront sensing. 

However we do not observe such effects when the 
optimal tomographic reconstruction matrix Wopt is used. 
Indeed, the analytical expression of this matrix at equation 
(|25p shows that it uses some a-priori about the turbulence 
profile as well as the noise, allowing to compensate for 
anisoplanatism and noise. 

One of the main justification of GLAO is to have 
reduced-seemg limited images. Although providing better 



angular resolution, the corrected PSF is still seeing-limited, 
meaning that it is not sharpened enough to gather the light 
into its inner regions and increase the spectroscopic SNR. 
As a conclusion, we have shown here that even if GLAO 
seems to be very promising, it is however not sufficient 
to allow to perform the 3D sp ectroscopy of distant galax- 
ies. iLe Louarn fc Hubinl ([200J) have suggested to combine 
GLAO with LGS, and found that using 4 LGS allow to im- 
prove the ensquared energy over a field of view equal to 
1 arcmin. However, this compensated field of view is still 
one order of magnitude lower that the 10 x 10 arcmin^ re- 
quired. This means that using GLAO with LGS over the 
Nasmyth FoV would again require at least one LGS per sci- 
entific target, and we have seen that such an architecture 
dramatically increase the cost of such an instrument. As a 
conclusion, GLAO is not usable in a system like FALCON, 
and the combination of MOAO with atmospheric tomogra- 
phy is the best architecture to reach our scientific goals. 

6.2 Influence of other parameters 

The goal of the results shown in the section [S] was to allow 
us to give some first order of specifications for the required 
FALCON AO system, in terms of sky coverage, limiting 
magnitude and number of actuators. We therefore found 
that using stars with magnitudes i? ^ 17 would allow us 
to satisfy our specifications with a sky coverage of at least 
50% for any galactic latitude. However, as we focused only 
on the spatial aspect of tomographic reconstruction, we 
did several simplifications before running our simulations, 
whose consequences are far from being negligible on the 
AO system's definitive design. 

First we assumed that we did not have any focal plane 
obstruction problem, meaning that the physical size of the 
AO components (the WFS or the IFUs) would not be a 
problem. For a real system, it is possible that the size of the 
optomechanical devices will be too large, implying that the 
first closest NGS to the scientific target will not be usable 
to perform wavefront sensing, and that it will be required 
to use more distant NGS. In fact, the equation (|28|) gives us 
the expression of the cdf Pi{D\ ^ r) of Di, the distance of 
the first closest NGS, and its expression is P\{D\ ^ r) = 
1 — exp— TTCT.r^. If we assume that all the devices have the 
same diameter Db, and that they can work even when they 
are tangent, then the distance D\ will indeed be the device's 
diameter Db, and the function Pi(Di r) will give us 
the probability Pn of not being able to use the first closest 
NGS. We therefore found that Pjv and Db are linked by the 
following relation: 

a, being a number of stars per solid angle, the 
above expression has to be multiplied by the focal plane 
scale (35 mm arcmin~^ for the VLT Nasmyth focus) to 
give a physical size. As an example, for b ~ —30°, we 
had a, — 0.35 stars per square arcmin {Rmax = 16) or 
(T» — 0.57 stars per square arcmin {Rmax ~ 17). Then, 
if we accept to loose 50% of the first closest NGS, we 
found respective diameters Db of 28 mm or 22 mm. These 
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Figure 17. Evolution of the median ensquared energy in a 0.25 X 0.25 arcsec^ square aperture for each galactic latitude b and each 
limiting magnitude Rmax in the J and H band, as a function of the number of corrected Zernike polynomials, in the case of GLAO 
correction (averaging of off-axis measurements). The full line corresponds to the median performance at the galactic pole, the dashed 
one to the median performance at b f» —60°, and the dash-dotted one at b « —30°. The dashed horizontal line shows the specification, 
i.e. 30% in J band and 40% in H band. 



numbers will of course be lower if we want to loose less stars. 

Another critical issue for AO is the temporal error due 
to time delay, which we did not take into account for our 
simulations, and which is going to degrade the performance 
of the AO system, in particular for high wind velocities. 
Additional simulations taking into account time delay are 
required to fix the t emporal bandw i dth o f the system. 
For MCAO systems, iLe Roux et al.l (|2004h have shown 
that the use of a Kalman filter allows to reach very good 
performance. Such a system therefore needs to be simulated 



to see if it could be incorporated into the FALCON's AO 
system. 

We have also to emphasize the fact that we simulated 
an open-loop system, where the off-axis wavefront sensors 
do not get any feedback from the on-axis DM. Moreover 
we assumed in our simulations perfect components with 
no errors and with instantaneous response. As we work 
in open-loop, the presence of errors like WFS aliasing 
or misalignents, DM hysterisis, and more generally non 
linearity, is going to have some non-negligible consequences 
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Figure 18. Evolution of the median full width at half maximum (in arcsec) for each galactic galactic latitude b and each limiting 
magnitude Rmax in the J and H band, as a function of the number of corrected Zernike polynomials, in the case of GLAO correction 
(averaging of off-axis measurements). The full line corresponds to the median performance at the galactic pole, the dashed one to 
the median performance at 6 Ri —60°, and the dash-dotted one at 6 Ri —30°. The dashed horizontal line shows the specification, i.e. 
FWHM = 0.25 arcsec. 



on the performance of the AO system. A tolerance analysis 

is therefore required to quantify the maximum errors we 
can accept for each AO subsystem. 

At last, the simulations we performed assumed static 
turbulence conditions, these latter being the median ones at 
the Cerro Paranal. However the turbulence is going to sig- 
nificantly change during the long exposures times required 
by our science objectives. Additional simulations with non- 
static turbulence conditions are therefore required in order 
to have some definitive realistic instrument performance. 



7 CONCLUSIONS 

We have proposed in this paper FALCON, a new concept 
of multi-object spectrograph for the ESO Very Large 
Telescope. Thanks to the combination of adaptive optics 
and atmospheric tomography methods, FALCON will allow 
to perform the 3D spectroscopy of several galaxies located 
up to 2 = 1.5 with a spatial resolution of 0.25 arcsec and 
a spectral resolution R = 10000, in a wide FoV. Such a 
performance implies to use Multi-Objects Adaptive Optics 
systems, and we have shown that for median Cerro Paranal 
atmospheric conditions, such systems will allow to reach 
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sky coverages of at least 50% for any galactic latitude. 

In terms of AO components, our simulations have 
shown that very sensitive wavefront sensors measuring 
the wavefront from guide stars with magnitudes _R ^ 17 
will be required, with 10 x 10 subapertures or actuators. 
The same number of actuators will be required for the 
deformable mirror in each integral field unit. Moreover each 
deformable mirror will require measurements from three 
oil-axis wavefront sensors, as atmospheric tomography 
methods are needed to improve the sky coverage. 

Such an instrumental concept is going to be a real 
technological challenge. Firstly because of the required 
miniaturisation of the AO components in order to be able 
to use in parallel several AO systems in the VLT focal 
plane and to avoid focal plane obstruction: the maximum 
size of the wavefront sensors and integral field units should 
not exceed 22 mm. Moreover, as atmospheric tomography 
is required in order to reach high sky coverages for any 
galactic latitude, each independent AO system will work in 
open-loop, meaning to use non-classical AO architectures. 
We have indeed shown that atmospheric tomography is 
definitely required for the required science objectives, as 
Ground Layer Adaptive Optics will not provide sufficient 
performance. 

Additional studies are required for the final design of 
the instrument, in particular for the temporal bandwidth of 
the AO system and the sensitivity to the variation of at- 
mospheric parameters. But the extrapolation of a concept 
like FALCON to the next generation of Extremely Large 
Telescopes is very promising. As an example, for a 42 me- 
ter telescope with a F/15 focal ratio, a 10 x 10 arcmin^ 
field will cover a physical size of w 2 x 2 square meters in 
the focal plane. A whole field instrumentation will therefore 
be impossible to manufacture because of tfie sizes of the 
optical components, and a MOAO system will be the best 
solution for wide field observations, which will be needed for 
the studies of the very early universe. Then, the use of an 
instrument like FALCON, but scaled to an ELT, will allow 
to detail galaxy physics with scales of 400 pc up to a = 7, 
and to understand how galaxies formed since the reioniza- 
tion epoch. 
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